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ABSTRACT 


Experimental data on radar response to soil moisture is analyzed with 
respect to the influence of soil texture. A truck-mounted scatterometer 
operating between 1 and 8 GHz, at incidence angles of 10° to 20° and at 
all linear polarizations observed test plots of sandy loam, silty clay 
loam and silty clay, Surface soil moisture conditions varied from the 
wilting point to field saturation. All test plots were non-vegetated 
and smooth with RMS surface variations * 1 cm. 

For a single soil texture the linear correlation between observed 
radar backscatter coefficient (dB) and measured gravimetric or volumetric 
soil moisture in the upper 5 cm of soil is typically 0.7 to 0.9. However, 
the sensitivity of radar response to moisture was found to be dependent 
upon soil texture such that sensitivity is inversely proportional to clay 
content of the soil. Thus, radar moisture estimation algorithms using 
either gravimetric or volumetric measures of soil moisture are expected 
to exhibit poor performance for complex multi -texture terrain if soil 
texture is unknown. Empirical evaluation of an equally distributed 
three texture ground scene produced linear correlations of less than 
0.7 for the upper 1 cm of soil and correlation was found to degrade drama- 
tically when considering larger depth intervals of soil (r:0.0 for the 
0-15 cm soil layer). 

In sharp contrast, estimation algorithms incorporating some knowledge 
of the tension at which soil water is retained by the soil produced strong 
correlation with sensor response, typically r = 0.8, regardless of depth 
interval considered within the surface 15 cm. The tested tension depen- 
dent soil moisture indicators include percent of field capacity, moisture 
normalized at soil tensions between 0.04 and 31 bars, volumetric moisture 
in excess of a transition moisture (as defined from dielectric investi- 
gations) and estimated soil tension. Importantly, these algorithms were 
found to be relatively independent of soil texture. 

The soil moisture indicators which produce radar estimation algor- 
ithms independent of soil texture are evaluated in terms of value to 
the potential user community of remotely-sensed moisture data. Normalized 
moisture indices will be difficult to interpret. Moisture in excess of 
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a transition moisture may prove to be of value if the transition moisture 
can be uniquely defined. Analysis of the radar data indicates that the 
transition moisture corresponds to an adsorbed w^ter layer at least 4-6 
molecular layers thick surrounding soil particles. Soil tension is found 
to be the most meaningful moisture indicator since it defines many critical 
points in the hydrodynamic behavior and agronomic importance of soil water. 
The current use of soil tension is mainly limited by the difficulty of 
measuring this parameter over the wide range of values typically assumed 
under field conditions. The high correlation between radar response and 
soil tension, a parameter not readily available even as a point measure- 
ment, should be carefully considered in any cost-benefit analysis of radar 
remote sensing of soil moisture. 


Key words: radar, scatterometer, soil moisture, soil texture, soil 

tension, remote sensing 
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1.0 INTRODUCTION 


In recent years, much interest has been expressed in developing 
capability for the operational remote sensing of soil moisture. Soil 
moisture information is needed on a local, regional and global scale for 
input into models in the fields of hydrology, meteorology and agriculture. 

In situ, point measurements of soil moisture are economically impractical 
in accounting for the temporal and spatial variability of soil moisture 
over large areas. Remote sensing research to date has investigated the 
correlation between soil moisture and surface albedo (Idso, et. al., 1975a), 
thermal inertia of terrain (Idso, et al., 1975b; Schmugge, et al., 1978; 

Idso, et al., 1976), radar backscatter coefficient (Cihlar and Ulaby, 1975; 
Batlivala and Ulaby, 1977) and the microwave emissivity of soils (Ulaby, 
et al., 1974; Newton, 1977; Schmugge, 1978). 

Hydrologists require knowledge of soil moisture in the upper 10 cm 
for the partitioning of rainfall into runoff, deep percolation and ground 
storage for assessment of watershed yield, reservoir management models 
and flood forecasting (Lawson, 1977; Salomonson, et al., 1975; Armstrong, 
1978). Regional meteorological models require soil moisture data for 
approximately the same surface layer to estimate moisture and heat fluxes 
into the atmosphere from land surfaces (Miller and Thompson, 1970). 
Agricultural interest in soil moisture ranges from that present in the surface 
layer which determines: tillage calendars, the germination and early growth 

of plants, the subsequent incidence of crop diseases, and hatching rates 
of insect populations (Idso, et al., 1975c; Hammond, 1975) to that present 
at depths of one or more meters in the soil which effects the growth of 
maturing crops and their yield (Erickson, 1977). At present, all remote 
sensing techniques are sensitive to the moisture present only in the surface 
layer, therefore moisture information at depths of one or more meters must 
be derived by modeling techniques (Baier and Robertson, 1966 and 1968; 
Hildreth, 1978; Reed, 1977, Cihlar and Ulaby, 1977; Rango, 1978). 

Events of the past decade have emphasized the fact that we live on 
a planet of finite resources and capability to efficiently regenerate 
renewable resources, The politics of food have important potential 
implications due to trends in world population and the vagaries of world 


climate and food production ('Schneider, 1976; CIA* 1976; Bryson, 1977; Smith, 
1976). Rational resource planning and development requires inventories 
and periodic monitoring of pertinent resource parameters such as soil 
moisture. 

Microwave sensors offer one of the most promising techniques for the 
remote sensing of soil moisture because of their all weather, day or 
night capabilities combined with an ability to penetrate a vegetation canopy 
with minimal attenuation of the emitted or backscattered signal. The great 
potential of microwave remote sensing systems is due to the large difference 
in the dielectric properties of water and those of dry soils at microwave 
frequencies. 

Recent research efforts have focussed primarily upon the selection of 
sensor parameters: frequency, incidence angle and polarization in a 

configuration that will minimize the impact of target parameters other 
than soil moisture. Pertinent target parameters include: vegetation type, 
density and moisture; surface aspect or slope; surface roughness; and soil 
type. Numerous experimental investigations using spaceborne, airborne 
and truck-mounted passive and active microwave sensors have demonstrated 
the feasibility of sensing surface soil moisture at microwave frequencies 
(Cihlar and Ulaby, 1975; Batlivala and Ulaby, 1977; Newton, 1977; Eagleman, 
1975; and others) and will not be discussed herein. Site specific experi- 
ments have verified the following as regards radar response to soil moisture 
(Batlivala and Ulaby, 1977; Ulaby, et al., 1977; Ulaby, et al., 1978; Ulaby 
et al., 1979; Bush and Ulaby, 1977): 

a) the radar backscattering coefficient is sensitive to soil moisture 
content at all microwave frequencies, 

b) signal attenuation by crop canopies is minimized at frequencies 
below 8 GHz, 

c) effects of surface roughness typical of dryland farming practices 
are minimized by the selection of a 4.5 GHz like polarized system operating 
between 7° - 17° incidence angle, 

d) sensitivity of radar response to surface slope is reduced at 
frequencies greater than 3 GHz, and 

e) radar sensitivity to soil moisture is dependent upon some aspect 
of soil type. 
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Soils are derived from a host of parent materials and subsequently 
modified by the interaction of complex geographic distributions of climatic 
regimes and associated plant communities. Soil associations vary regionally 
as a function of the above parameters, and soil texture as defined by the 
particle size distribution of soil solids is observed to vary considerably 
over tens of meters. Thus, assuming that other target "noise" such as 
water bodies and man-made structures can be filtered from the radar data 
and that vegetation and roughness effects are minimized by optimizing 
sensor selection, effective soil moisture estimation algorithms must account 
for the observed dependence of radar response to soil type. 

1.1 Statement of Problem and Purpose 

The complex dielectric constant and thus passive and active sensors 
operating at microwave frequencies have all been observed to exhibit some 
dependence upon soil texture in their response to soil moisture. Since 
most prior experimentation was site specific and examined only one soil 
texture, the influence of soil texture upon the sensor response to soil 
moisture was found to be significant only through comparison of results from 
several experiments. 

The purpose of this study is to: a) experimentally verify prior 
observations that 1-8 GHz radar response to soil moisture is influenced by 
soil texture, b) investigate the behavior of soil water as a function of 
soil texture and the electromagnetic characteristics of the soil solution 
at microwave frequencies, c) examine available quantifiers of soil moisture 
such as gravimetric and volumetric moisture, percent of field capacity and 
soil tension with respect to radar response and soil texture, and d) evaluate 
those moisture indicators found to produce texturally independent soil 
moisture algorithms in terms of their value to a user community of soil 
moisture information. 

In order to accomplish the above objectives, the 1977 Bare Soil Experi- 
ment (Dobson, 1979) was conducted using the University of Kansas truck- 
mounted Microwave Active Spectrometer 1-8 GHz system to observe carefully 
selected and prepared bare soil fields of sandy loam, silty clay loam and 
silty clay at incidence angles between 10° and 20° relative to nadir, 
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Radar backscattering measurements were made at each field for surface 
soil moisture conditions ranging from the wilting point of plants to 
saturation and partial flooding of the test fields. In conjunction with 
the radar measurements, many target parameters were monitored during the 
course of the experiment including: gravimetric moisture profile, soil 
bulk density, soil texture, organic matter content, surface slope, surface 
roughness, soil water desorption characteristics and precipitation. 

The radar data acquired from this experiment was used to develop 
algorithms using various soil moisture indicators for each soil texture 
independently and for all soil textures considered simultaneously. These 
algorithms were then compared to evaluate their respective independence to 
soil texture. 

1.2 Chapter Outline 

This study is divided into five chapters. Chapter 2 treats soil 
textural effects upon soil water retention and upon the electromagnetic 
behavior of soil water at microwave frequencies. Soil water retention is 
examined with respect to soil composition, structure and mineralogy which 
collectively determine the pore size distribution, the specific surface 
area and the affinity of the soil solution to the soil solid surfaces. 

The behavior, composition and structure of the soil solution is examined 
as a function of distance from soil particle surfaces. Dielectric measure- 
ments of moist soils at microwave frequencies are reviewed and linked to 
results of soil moisture experiments utilizing radiometers and radar for 
observation of a variety of soil textures. Soil moisture quantifiers are 
reviewed with respect to their dependence on soil texture. 

Chapter 3 describes, in brief, the 1977 Bare Soil Experiment. Target 
characteristics and the homogeneity of each test field are summarized. The 
angular and frequency response of the radar with respect to gravimetric 
soil moisture is examined and found to be consistent with results derived 
from other experiments with the MAS 1-8 GHz system (Batlivala and Ulaby, 
1977). 

The results of the 1977 Bare Soil Experiment are analyzed in Chapter 
4 to evaluate the textural independence of the following soil moisture 
indicators with respect to radar backscatter: gravimetric moisture, 
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volumetric moisture, normalized moisture (such as percent of field capacity), 
volumetric moisture above some transition moisture and soil tension. 

Chapter 5 presents an evaluation of those soil moisture indicators 
found to show minimal dependence upon soil texture in the preceding 
analysis which can be readily Interpreted by the potential user community 
of agronomists, hydrologists and meteorologists. Specific recommendations 
are made for future experiments in light of these findings. 

2.0 BASIC CONSIDERATIONS OF SOIL WATER AMD THE ELECTRICAL BEHAVIOR 

OF SOIL WATER 

As a prelude to an understanding of the interaction of an impressed 
microwave field on the soil -water system, it is necessary to address the 
physics of the soil and its components. Soil consists of unconsolidated 
mineral material derived as the weathering product of parent rock either 
in situ or as a sediment of some transporting agent (water, wind or ice). 

Soil can be considered as a three component system: soil solids, soil 
water and air; the interaction of these components imparts characteristic 
structure, soil water and density relationships within a given soil. 

In the surface horizon of agricultural soils, soil solids consist 
of organic and inorganic mineral fractions. The organic fraction is 
concentrated near the surface and consists of live plant material, such 
as roots, decaying organic debris and carbonaceous decay resistant humus. 

The organic fraction is usually less than 10% of the dry soil by weight 
(and typically less than 5%); it is extremely active in terms of soil 
chemistry and soil moisture dynamics, especially the colloidal humus. 

Soil organic matter is a prime determinant of surface soil structure, 
the water-holding capacity of soils and the soil's cation exchange 
capacity. The inorganic fraction consists of the weathered mineral 
grains of the parent rock. The grains can be described in terms of 
their mineralogy and particle size distribution. 

2.1 Soil Moisture Concepts 

Soil chemistry and water-capacity are governed primarily by the 
interaction of soil structure and soil composition. Soil structure 
establishes the volume and size distribution of pore spaces which can 
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be occupied by soil water. Soil composition, defined herein as the 
soil particle size distribution and particle mineralogy, determines the 
net effective soil surface area to which water molecules may adhere 
and the electro-chemical properties of the particle surfaces. 

2.1.1 Soil particle size 

The arrangement of solid particles determines the structure, norosity 
and bulk density of the soil. Soil solids range in size from colloids 
on the order of 50 A 0 in diameter to gravel. Soil particle sizes are 
classified into categories of soil textural separates. Many classification 
schemes are available (Figure 2,1), but use of the 1951 United States 
Department of Agriculture (USDA) size categories will be adhered to in 
the present discussion. Soil particles are classified into textural 
separates by mechanically passing the inorganic fraction (either wet 
or dry) through sieves and subsequent sedimentation of the fine material. 
Sieving and sedimentation techniques both assume spherical soil particles; 
however, soil particles are not generally spherical and this assumption 
is particularly weak for the clay size fraction. In mechanical sieving, 
the mesh size of the screen refers to the equivalent diameter of a sphere 
that will pass through that screen. Textural separates smaller than 
sand must be isolated using sedimentation techniques based upon Stoke's 
law which predicts the settling velocity of spherical particles for a 
given set of environmental conditions. 

V » 2/9 gr 2 fc . ^ -) (2.1) 

where: V = the steady falling velocity of a particle 

g = acceleration due to gravity 
r = particle radius 
p s = particle deni sty 
p w = liquid density 
z = viscosity of liquid 

Soil particles are classified as having diameters equivalent to 
those of spheres of the same denisty and settling velocity. The sedimentation 
technique commonly employs either a hydrometer to monitor the density of 
the soil suspension as a function of time or pipetting of the suspension 
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2.1 Classification of soil separates $ 2.0 mm on the basis of particle si 2 e (from 
Baver, Gardner and Gardner, 1972). 



to gravimetrically determine the mass of soil solids remaining in 
suspension as a function of time. It should be remembered that all 
determinations of soil textural separates are biased by the relative 
validity of assumed spherical particles. The sphericity of a particle 
is governed by its mineralogy and weathering history and will be more 
fully discussed in the following section. 

Since all so/ls contain a distribution of particle sizes, it is 
convenient to classify a soil by the weight percent of the soil with- 
in specific size categories (sand, silt, clay). For a soil having the 
particle size distribution shown in Figure 2.2, the soil textural separates 
of sand, silt and clay represent the integrals of the distribution 
between the limits shown ratioed to the integral of the entire distribution. 
The percent by weight of each size interval then determines the soil 
textural class using the USDA triaxial classification based on sand, silt 
and clay fractions shown in Figure 2.3. 

2.1.2 Soil Particle Surface Area 

The net soil surface area available to interaction with the soil 
solution consists of the sum of component particle external surface 
areas and clay particle internal (interlayer) areas as limited by soil 
particle contacts. 


S = 2A e + sA i (2.2) 

2 

where;; S = the specific surface of a soil in m /g 

2 

A g - external surface area of a soil particle in m /g 

2 

A. = internal surface area of a soil particle in m /g 

The external surface area of soil particles is a function of particle 
size and shape. Figure 2.4 shows the relationship between particle size, 
the number of particles per gram and external surface area assuming 
spherical partirles. It is apparent that .002mm diameter (2 p) clay has 
approximately 50 times the external surface area of an equal weight of 
.1 mm diameter fine sand. Table 2.1 assumes closest hexagonal packing of 
spherical particles to compare external surface area on a volume basis. 

It is obvious that, for a given soil, A g is dominated by the percent clay 
fraction. 
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TABLE 2.1 


The relation of particle size to external surface 
(from Baver, Gardner and Gardner, 1972). 


Diameter pf sphere 

Textural name 

Volume per particle 
(n/6) D 2 

Number of particles 
m > 6 - cc 

Total surface 
nD 2 x number of particles 

1 cm 

Gravel 

(tr/6) (1)3 

l 

3.14 

cm 2 (0.49 in. 2 ) 

0.1 cm 
{. mm) 

Coarse sand 

(n/6) (O.l) 3 

1 x 10 3 

31.42 

cm 2 (4.87 In, 2 ) 

0.05 cm 
(500 V) 

Medium aand 

(n/6) (0.05)3 

8 x 10 3 

62.83 

cm 2 (9.74 in. 2 ) 

0.01 col 
(100 ») 

Very fine sand 

(a/6) (0.01)3 

1 X 10 6 

314.16 

cm 2 (48.67 in. 2 ) 

0.005 cn 
(50 u> 

Coarsu Silt 

(tr/6) (0.005) 3 

8 x 10 6 

628.32 

cm 2 (97.34 in. 2 ) 

0.002 ca 
(20 |i) 

Silt 

(n/6) (0.002)3 

125 X 10 6 

1,570.8 

cm 2 (1.69 ft 2 ) 

0.0005 cn 
(5 h) 

Pine silt 

(n/6) (0.0005)3 

8 x 10® 

6,283.2 

cm 2 (6.76 ft 2 ) 

0.0002 cm 
(2 U) 

Clay 

(n/6) (0.0002)3 

125 X 10® 

15,708 

cm2 (16.9 ft 2 ) 

0.0001 cn 

(1 u) 

Clay 

(tt/ 6) (O.OOOl) 3 

1 x IQ 32 

.31,416 

cm 2 (33.8 ft 2 ) 

0,00005 cn 
(500 mu) 

Clay 

(n/6) (0.00005) 3 

8 x 10 1? - 

62,832 

cm 2 (67,6 ft 2 ) 

0.00002 cm 
(200 mu) 

Colloidal Clay 

(n/6) (0.00002)3 

125 x 10 12 

157,080 

cm 2 (169 ft 2 ) 

0.00001 cm 
(100 mu) 

Colloidal Clay 

(n/6) (0.00001) 3 

1 x 1015 

314,160 

cm 2 (338 £t2> 

0.000005 cn 
(50 mU) 

Colloidal Clay 

(n/6) (0.000005) 3 

8 x 10 15 

62B.320 

cm 2 (676 ft 2 ) 
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Many soil characteristics are surface phenomena: particle attraction, 

water, gas and nutrient adsorption among others. Figure 2.5 emphasizes 
the significance of the high specific surface area of clay sized 
particles. 

The rate of increase in specific surface area with decreasing particle 
diameter is amplified by considerations of particle shape. The spherical 
particles assumed above have the smallest possible area to volume ratio; 
as particles become less spherical, the area to volume ratio increases. 

Soil particles are rarely spherical and, in general, become more plate- 
like as effective diameter decreases. The loss of sphericity with 
decreasing particle size is the result of mineralogic control as 
influenced by weathering processes. Sand and most silt sized particles 
are composed of primary weathering minerals, usually quartz, feldspars 
and ferro-magnesium minerals which may have been mechanically reduced 
in size; these particles are subangular to rounded in shape depending, 
in part, upon abrasion. Silt size particles are diverse, irregularly 
shaped fragments which are rarely smooth or flat; similar to sand, 
quartz is the dominant mineral. Clay sized particles include primarily 
clay minerals mixed with relatively coarse grained amorphous crystalline 
minerals (Table 2.2). The clay minerals are complex layer-lattice 
alumi no-silicate sheet structures which have a characteristic plate-like 
shape and are often colloidal. Soil clay platelets may be only 50 A 0 
in diameter and 10 A 0 thick where lp = 10^ A 0 (Taylor and Ashcroft, 

1972). 

Table 2.3 demonstrates the effect of particle shape on external 
surface area. The introduction of plate-like particle shape produces 
a dramatic increase in external surface area when compared to a 
reference sphere. A disk .2p thick has a surface area 257 times greater 
than the 2p reference sphere where 2p is the maximum equivalent diameter 
of clay size particles. 

Internal surface area, A • , is primarily a phenomenon of the clay 
minerals. The clay minerals can have very high internal surface areas 
in the interlayers of the particle crystal lattice. It is necessary to 
introduce some basic concepts of clay mineral structure in order to 
understand the nature of clay internal surface area. 
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TABLE 2,2 

Principal Secondary Minerals of the Clay Size Fraction 
(from Townsend , 1973). 



oxides 

silica SiOa.nHaO 

hydrated iron oxides Fe(OH)i.nHsO 
alumina AhOa.nHaO 

Amorphous 

silicates 

allophano AIjOj . 2SiOa . nHaO 
hisengerito FeaO»,2SiOs,nHjO 

phosphates 

evansite Al3PO-»(OH)o,nHsO 
azovskite FesPO^OHJa.nHaO 

* 

carbonates 

calcite CaCOs 
magnesite MgCOa 
dolomite CaCOa. MgCOa 
siderite Fe COa 

Crystalline 

oxides 

goethite aFeO.OH 
lepidocrocite yFeO.OH 
maghemite rFeaOa 
gibbsite AkOa 


silicates 

complex alumino-silicates, grouped 
according' to their structural 
characteristics. 












TABLE 2.3 

Surface Area -In Relation to Shape of Particle 
(from Baver, Gardner and Gardner, 1972). 


Stull’ i* 


Radius 

(cm) 

Vnlur.e 

(ca?) 

Surface 

Increase 
in surface 
(percent:) 

Sphere 


1 X 10-'* 

4.2 x 10~ 12 

1.26 x 10 -7 


Disk , 

h - 1 x 1(P 

cm 

1.155 x 10" 4 

4.2 x 10 -12 

1.56 x 10~ 2 

23. B 

h - 5 x 10 -5 

cm 

1.67 x 10 -4 

4.2 x 10~ 12 

1.84 x 10" 7 

45.8 

h « 2 x 10" 5 

cm 

2.58 x 10* 4 

4.2 x 10“ l2 

4.51 x 10’ 7 

257.8 

h - 1 x 10~ 5 

cm 

3.65 x 10" 4 

4,2 x 10 -12 

8.59 x 10 -7 

538.9 
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2.1.3 Some Properties of Clay Minerals 

Since quartz particles are very Inactive chemically, the sand and 
silt sized soil fractions contribute little to soil chemistry. On the 
other hand, the clay minerals are quite surface active due to their 
aluminosilicate sheet structure and dominate the chemical behavior of 
the soil inorganic fraction. 

Based upon mineral structure, three main types of aluminosilicate 
layer clays can be identified: montmoril Ionite, illite and kaolinite. 

Table 2.4 lists other principal vystalline silicate clay minerals as 
well. The distinguishing features between various clays are the 
numbers and kinds of layers in a crystal unit layer and the isomorphous 
substitutions of ions within these structures. 

The basic molecular components of clay aluminosilicates are the 
silica tetrahedron and the alumina octahedron, a silicon atom surrounded 
by four oxygen atoms and an aluminum atom surrounded by six hydroxyls or 
oxygen atoms respectively. Each basic building block or crystal unit cell 
is arranged in planar sheet structures bound by common oxygen atoms. The 
tetrahedral and octahedral sheets are variously arranged and stacked, and 
are bound by common oxygen atoms to produce the crystal unit layers dis- 
tinctive of each aluminosilicate layer clay. Thus, crystal unit layers 
of the 1:1 family consists of one silica tetrahedral sheet bound to one 
alumina octahedral sheet and the 2:1 family has an alumina octahedral layer 
in between two silica tetrahedral layers. The crystal unit layers are then 

stacked with various bonding energies in either single mineral form or 
as mixed layer clays to form the clay particles of an equivalent diameter 
as identified by mechanical size fractionation. 

The aluminosilicate layer clays are further classified by the 
isomorphous substitution of ions within the lattice of a given crystal 
unit cell and the bonding characteristics between unit layers. The 
idealized clay structures presented above are electrostatically neutral, 
this rarely occurs in nature. The idealized crystal unit layers are 
disrupted by the isomorphous substitution of other cations of similar 
size and lower valence in the crystal lattice which causes the crystal 
to become negatively charged. Table 2.5 gives the ionic radii of elements 
commonly substituted and adsorbed to silicate clays. These negative 
charges are commonly balanced by the presence of ah additional cation 
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TABLE 2.4 

Principal Crystalline Silicate Clay Minerals 
(from Townsend, 1973). 

I 



TABLE 2-5 

Ionic Radii* of Elements Common in 
Silicate Clays and an Indication of 
Which are Found in the Tetrahedral 
and Octahedral Layers (from Brady, 


Note that aluminum can fit in either layer. 


Ion 

R» Jli!» (A)* 

Found m 

Al»* 

0 411 
0.50/ 

Silica t.HMhfda 

Fe‘* 

064 


Mg** 

Zb** 

065 

070 

Alumina oclahcilra 

Fe** 

0.75) . 


Ca** 

Na* 

0.94 ( 
0,98 f 

Eacbange *ii« 

K* 

I.VJ 


O " 

1,45 

Both layers 
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adsorbed to the surface of the crystal layer in the vicinity of the 
crystal lattice substitution. The adsorbed cations are termed ex- 
changeable cations since they can be replaced at the surface of the 
clay layers by other cations and water molecules. The extent and 
nature of isomorphous substitution determines the net surface charge 
of the layer silicate clays. 

Kaolinite is the most prominent member of the 1:1 lattice clays. 

The alumina octahedral and silica tetrahedral layers are tightly held 
together by hydrogen bonding of pairs of adjacent oxygen and hydroxyl 
groups. The hydrogen bonding gives the crystal structure great strength 
and permits little substitution within the lattice; therefore cations 
and water are not absorbed between crystal units and no swelling of this 
clay ordinarily occurs. The specific surface of kaolinite is thus 
defined by its external surface area only, which is relatively small 
since the hydrogen bonding results in large particle sizes, typically 
.2 to 2 pm across (Brady, 1974). Hence, the surface activity of kaolinite 
as defined by the area density of exchangeable cations is low since 
substitution is rare and surface activity is primarily determined by 
broken bonds at the edges of crystal sheets in the lateral dimension 
along the particle edges. 

Hydrous micas are 2:1 lattice clays which, like kaolinite, are 
also non-expanding clays. Illite is the most common member of this 
group and is often found associated with montmoril Ionite. The clay 
mineral fractions of the soils examined during the 1977 Bare Soil Experi- 
ment are approximately 30% illite and 60% montmoril Ionite. Isomorphous 
substitution occurs mainly in the silica tetrahedral layer where approx- 
imately 15 to 25% of the tetravalent silicon atoms are replaced by larger 
trivalent aluminum atoms (Townsend, 1973). Non-exchangeable potassium 
ions on the surfaces of the crystal unit cell are tightly bound between 
adjacent tetrahedral layers, thus bonding adjacent unit cells. The 
potassium acts to restrict the expansion potential of the illite particle 
structure thereby limiting cation exchange sites to the particle external 
surface area. Cation exchange can occur at the planar surfaces of the 
particle and along the frayed lateral edges of the structure. Illite 
particles tend to be intermediate in size between those of kaolinite 
and montmorillonite. 
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Montmorll Ionite Is a 2:1 lattice clay with considerable potential 
to expand, Isomorphous substitution occurs in both the tetrahedral 
and octahedral layers. Commonly, aluminum substitutes for silicon in 
the tetrahedral layers and magnesium and iron substitute for aluminum 
in the octahedral layers. The high degree of substitution imparts a 
large negative net surface charge on montmorillonite crystals which is 
satisfied by a swarm of exchangeable cations adsorbed to both the 
internal and external surfaces of a particle. Within a particle, the 
clay crystal unit cells are very loosely held by only van de Waal's 
forces and weak oxygen to oxygen affinities, This weak bonding of the 
clay unit cells allows the tremendous expansion of the particle 
structure when cations and water are introduced into the interlayers of 
the clay lattice. The very large internal surface of montmorillonite 
combined with a large relative external surface area (typically .01 
to 1 \m equivalent particle diameter) and the high degree of 
isomorphous substitution causes the cation exchange capacity (CEC) of • 
montmorillonite to be an order of magnitude greater than that for 
either kaolinite or illite (Table 2.6). 

In summary, the total CEC of a given clay mineral is determined 
by the degree and type of isomorphous substitution within the crystal 
lattice and the number of broken bonds at the frayed crystal edges, 

The part of the CEC determined by isomorphous substitution is relatively 
constant in magnitude for a given soil, while the CEC determined by 
broken edge bonds is dependent upon the pH of the soil solution. At 
the crystal edges, the contribution to total CEC is proportional to 
soil pH; as pH increases more hydroxyl groups are dissociated with the 
hydrogen being readily .exchangeable. The net effect of this pH 
dependent charge is small compared to the particle surface charge 
determined by isomorphous substitution (Figure 2.6), 

The CEC and the relative specific surface of clay minerals is 
important to the present discussion because of the manner in which 
electrostatically charged solid surfaces modify the distribution and 
electromagnetic behavior of water molecules in the soil solution. Also, 
it is obvious from Figure 2.5 that the cation exchange capacity of even 
small quantities of the colloidal fraction of organic humus will have 
a profound effect. These issues will be discussed after consideration 
of soil v/ater. 
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TABLE 2.6 

Comparative Properties of Three Major Types of Clay 
{from Brady, 1974). 


Type ofCIsy 

Propciljr Mommorillonile lliilc Kaolmhe 


Sire 

Shape 

Specific surface (m'/s) 
Exlcrn.il surface 
internal surface 
Cohesion, plnsticily 
Swelling capacity 
Cation exchange canxcily 
(meq/IOOg) 


0.01-1,0 0, 1-2.0 0.1 -5.0 

Irregular Oakes Irregular Bakes I lengonat crystals 

700-800 I 00-120 5-20 

High Medium Low 

Very high Medium None 

High Medium Low 

High Medium Low 

80-100 15-40 3-15 


ORIGINAL 

OF POOR 


PAGE is 
QUALITY 
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Cation Exchange Capacity (meq / lOOg) 



Figure 2.6. Influence of pH on the cation exchange capacity (a measure 
of the negative charge) of montmoril ionite and humus. Note 
that below pH 6.0 the charge for the clay mineral is relatively 
constant. This charge is considered permanent and 'is due 
to ionic substitution in the cryst'l unit. Above pH 6.0 the 
charge on the mineral colloid increases because of ionization 
of hydrogen from exposed 0-H groups at crystal edges. In 
constrast to the clay, essentially all of the charges on the 
organic colloid are considered pH dependent (from Brady, 1974). 



2.1.4 Measurement of Soil Water 

Water content of soils may be determined by direct or indirect 
methods. Direct methods use a variety of processes to separate water 
from the soil solids, such as the oven drying of soil, for comparison of 
the mass of water to the mass of soil. Indirect methods measure a 
property of the soil which is related to water content such as: gjunma 
ray and neutron attenuation, electrical conductivity and resistivity 
or the tension at which soil water is held. The last three methods are 
related to moisture content by its relationship to soil moisture tension, 
which is a pressure or work per unit volume usually referenced to 
atmospheric pressure in bars. 

All determinations of soil water are made with respect to some 
arbitrary reference zero. In the gravimetric determination of soil 
moisture, the reference zero is commonly the mass of soil after drying 
at 100 to 110°C to a constant weight, usually requiring about 24 hours. 
Other reference zero points are possible since there is nothing unique 
about dehydration at 110°C (Figure 2.7). It is necessary to understand 
the structural nature of water retained by soils to explain the observed 
differences in dehydration as a function of temperature. 

The inclusion of organic matter in the soil matrix further confounds 
the selection of a reference zero. The organic fraction of soil contains 
many volatiles which oxidize at temperatures of 110°C although the net 
effect is usually quite small for agricultural soils where organic matter 
content is often less than 5% of the total dry weight of soil. After a 
10 day drying period, a silt loam with 3 % organic matter was found to 
have a weight loss of; only 0.3% moisture by weight (Black, 1964). 

The volumetric moisture content of the soil, M y , is the volume 
fraction of the soil occupied by water. M y is a function of the weight 
percent of water and the bulk density of the soil , 


M v - .01 Mg X P B /P K 


where: M„ 


w 


« percent of water per gram of dry soil 

3 

= bulk density of the soil, g/cm 

3 

= density of water =1.0 g/cm 


(2.3) 
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Weight, Arbitrary Units 



Temperature, °C 

Figure 2.7. Weight loss of clay minerals as a function of temperature. 

The weight of the sample at a temperature of 800°C or 
greater is taken to be unity. Water content percentages on 
a mass basis may be obtained by subtracting 1 from the 
ordinate figure and multiplying by 100. (The numbers to 
the right of the decimal thus may be regarded as percent- 
ages.) (From Black, 1964). 
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In clay soils, gravimetric moisture M g and bulk density are not 
independent due to the expansion of clay particles by adsorption of 
interlayer water. Thus, is inversely proportional to p b in clay soils. 
The effects of this dependence can produce a considerable error in the 
calculation of volumetric moisture content from M g , up to approximately 
20%, if bulk density is treated as a constant as a function of M g . As 
a consequence, accurate volumetric moisture determination necessitates 
the concurrent measurement of soil bulk density of clay soils. Since the 
field measurement of bulk density is difficult and the accuracy of a 
particular method is often restricted to a small range of M , it is often 
impractical to make concurrent determinations of M g and p b . Appendix A 
further treats the problem of bulk density sampling and presents a simple 
physical model for estimating p b as a function of M g for clay soils. 

Soil is composed of solid, liquid and gas phases. At most, the 
liquid phase can only occupy that volume fraction of the soil not 
occupied by solids, or the pore space. The pore space of a soil is 

defined By soil bulk density and the specific density of the soil solids. 

The bulk volume of soil is given as: 

V b - M s /p b (2.4) 

3 

Where: V b = bulk volume of the soil, cm 
M s = dry mass of soil, g 

3 

P b = bulk density of soil, g/cm 

Similarly, the particle volume is determined by: 

Vp - H s /P p (2.5) 

3 

where: V p = the total volume of solid particles, cm 

M g - dry mass of solid particles, g 

3 

P - the specific particle density of soil ~ 2.65 g/cm 

r 
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From these definitions, it follows that the volume of soil occupied by 
solids is given by the ratio V p/V b ; and the porosity of the soil is 
1-V /V The total percent porosity S t can then be defined as: 

s t = 100 (1 - p b /p p ) (2.6) 

By definition, percent volumetric moisture must be less than or equal to 
porosity, 

100 M y * S fc 

It should be noted that for expanding lattice clays, p b is known 
to be dependent upon soil water content (Figure 2.8), Thus it is to 
be expected that the total porosity of clay soils will be proportional 
to moisture content. 

While maximal water content of a soil is determined by its porosity, 
the water content of a soil, at any given point in time, is determined by 
the interaction of soil structure and soil -water energy relationships. 

Structurally, soil water consists of pore water, solvate water, 
adsorbed water and structural water (Kezdi, 1974). Pore water occupies 
the macropores of the soil matrix; its volume is established by the 
effect of hydrodynamic forces acting upon the capillarity of the pore 
size distribution within the soil. The behavior of pore water is physically 
and chemically equivalent to that of bulk water having the same salt 
concentration. Solvate water forms a hydration layer around soil particles. 
The thickness of the solvate water layer around an individual particle is 
dependent upon the electrostatic charge of the particle surface, ionic 
bonding forces within the soil solution and the polar nature of water; 
the solvate layer is typically 20 to 400 A° thick which corresponds to a 
water layer 8 to 160 molecules thick depending upon the specific particle 
and solution considered (Bolt and Bruggenwert, 1976). Solvate water may 
hbve a greater density and viscosity than bulk water, although much of it 
can still flow under the hydrodynamic forces governing the flow of pore 
water. Adsorbed water is a tightly held hydration shell of 1 to 10 
monomolecular layers of water adjacent to the solid surface depending 
upon the mineralogy of the particle and the ion concentration of the 


A, Suction (cm of Water) 



Figure 2.8. The moisture characteristic and swelling curve of a 
sample of kaolin clay (from Childs, 1969). 
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solution (van Olphen, 1963). The adsorbed water cannot be moved by the 
influence of normal hydrodynami c forces and cannot be completely removed 
by traditional oven drying procedures at 110°C. Structural water is not 
true water in the formal sense of the term, and is composed of the 
hydroxyl groups within the crystal lattice structure of primarily the 
layer aluminosilicates. Structural water cannot be removed from the 
crystal lattice without severe modification of the crystal. 

The dehydration curves of common soil minerals in Figure 2.7 
illustrate the relative differences in the volume of water held in given 
structural states. The pore water and most solvate water is evaporated 
from all the minerals at 100°C. The remaining solvate water and the 
adsorbed water are typically evaporated at temperatures less than 200°C. 
The amount of structural water is shown by the relatively horizontal 
portions of the curves at 16%, 6% and 1% for kaolinite, montmoril Ionite 
and quartz respectively. Evaporation of hydroxyl groups is responsible 
for the weight loss of the clay minerals at about 600°C. Differential 
thermal analysis of a calcium montmorillonite (Figure 2.9) illustrates 
the amount of heat required to liberate the various forms of water held 
by the clay. All pore and solvate water is evaporated by the adsorbed 
heat between 0° and the endothermic peak at s 120°C. The shoulder on 
this endothermic peak is due to the vaporization of adsorbed water at 
temperatures between s 150 to 200°C. The loss of hydroxyl water at 
600°C causes a second endothermic peak. The exothermic peak at 900°C 
is attributed to phase changes in the crystal lattice structure and is 
preceded by the endothermic evaporation of the final traces of structural 
hydroxyl water (van Olphen, 1963). 

2.1.5 Energy of Soil Water 

The free energy of soil water characterizes the summation of 
potential, kinetic and electrical energy acting upon the soil water 
system. The dominant forces affecting the free energy of soil water are 
gravitational, matric and osmotic. The status of free energy of soil 
water is measured as a potential, where soil water potential is defined 
as the difference between the energy of soil water and that of pure 
water (Brady, 1974). The soil water potential is commonly referenced 
to the energy of pure water at an arbitrary water table. The soil 
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Figure 2.9. 


Differential thermal analysis (DTA) of a montmoril Ionite (schematic) (from van Olphen, 
1963). 



°c 


Endothermic Reaction 1 = Loss of Surface Hydration Water 
Endothermic Reaction 2 = Loss of Lattice Hydroxyl Water 
Exothermic Reaction 3 = Crystal Phase Change 


solution may occur as thin films covering the soil particles or as pore 
water at higher moisture contents. Only at moistures near saturation 
does the soil solution exhibit the properties of flow, vaporization and 
internal energy that characterize bulk water. As a result, the energy 
relations of film and pore water have proven themselves more useful in 
the study of pi ant-soil -water relations than the properties of bulk 
water (Taylor and Ashcroft, 1972). 

Total soil water potential is the sum of the potentials arising 
from each of the forces acting upon the soil water system (Figure 2.10). 
The gravitational force produces a positive potential and functions to 
transport water to lower elevation in the soil profile. The gravitational 
force acts primarily on the pore and some solvate water in the capillary 
maze of the soil matrix. The matric potential measures the attraction 
of the soil water molecules to each other and to soil solid surfaces 
due to cohesive and adhesive forces respectively. The cohesive force 
produces the surface tension of water. The phenomenon of capillarity 
is due to the combined effects of surface tension and adhesion of water 
to the soil solid surfaces. Thus, the matric potential is negative since 
it functions to reduce the free energy of soil water. 

The osmotic potential of soil water is also negative; it is the 
result of attractive forces of ions and other solutes which tend to 
hydrate by orienting polar water molecules. The electrostatic fields 
of ions cause a polarization and binding of surrounding water molecules 
in a manner described by a model of three concentric water regions 
enveloping the ion (Kavanau, 1964). An innermost region of structured, 
polarized and electrostructured water molecules grades into an inter- 
mediate region where the water is less structured than ordinary water 
and an outermost region of bulk water. Polarization is found to be 
proportional to the charge of the ion and inversely proportional to 
ionic radius squared (Nielson, et al , , 1972). Thus, the polarizing 
effects are observed to be most pronounced for small highly charged 
ions such as sodium (Shainberg and Kemper, 1966). 

Since osmotic and matric potentials are negative, it is convenient 
to describe them in terms of positive tensions or suctions by which the 
water is held within the soil. The tension of soil water describes the 
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Figure 2,10. Soil water potential. In an isothermal equilibrium system, 
matric suction is the pressure difference across a mem- 
brane separating soil solution (equilibrium dialyzate) 
in bulk from soil containing the solution, the membrane 
being permeable to solution, but not to matrix or to 
mass flow of gas; osmotic suction is the pressure differ- 
ence across a semi permeable membrane separating bulk 
phases of pure water and the soil solution; total suction 
is the sum of the matric suction and theosmotic suction 
and is the pressure difference across a semipermeable 
membrane separating pure water and soil that contains 
solution. An ideal semipermeable membrane is permeable 
to water only (modified from Black, 1964). 
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free energy of the soil water and its availability to plants or for 
mass movement under the influence of hydrodynamic forces. Suction is 
commonly expressed in terms of the height of a water column whose weight 
equals the suction or as pF which is the logarithm to the base 10 of 
this column height in cm. Approximate equivalents between suction 
expressed as cm of water, pF and atmospheres are given in Table 2,7. 

In the American literature, suction or tension is more commonly expressed 
in terms of atmospheres or bars of pressure. The relevant conversion 
factors are: 

fi ? 

1 bar - 1 x 10 dynes/cm 

1 millibar =1 x 10 3 dynes/cm 2 

1 bar = 0.98692 atmosphere 

1.0133 bar ** 1.0 atmosphere 

1 bar = 1033.26 g/cm 2 

1 bar « 1022.7 cm of pure water at 25°C 

Soil water energy potential, as measured by soil suction, is related 
to soil water content for a given soil and moisture history. However, 
soil tension is not a direct expression of the thickness of water films 
around soil particles but also depends upon soil structure, specific 
surface and the solute content of the soil water. 

The differences in soil water retention at various tensions between 
soil textural classes observed in Figure 2.11 correspond to differences 
in pore size and shape distributions combined with differences in the 
nature and strength of adsorptive forces for soil particles with differing 
mineralogy. Coarse textured soils, such as sand, are characterized by 
comparatively large pores which can hold water only at low suctions. On 
the other hand, clay soils have a much finer network of pores and soil 
water retention is dominated by the interaction of soil water with a 
large relative surface area. Thus, a considerable range of the -water 
content of clays is highly controlled by higher tension matric and osmotic 
forces . 

An explanation of the dramatic differences in soil water retention 
between soil textural classes at suctions greater than ~ 1 bar lies 
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TABLE 2.7 

Approximate Equivalents of Common Means of Expressing 
Differences in Energy Levels of Soil Water 
(from Brady, 1974). 


Height of Unit Column 
of Witter 
(cm) 

Logarithm of 
Water Height 
(PF) 

Atmosphere (bars) , 

10 

1 

001 

100 

2 

0.10 

346 

2,53 

1 

i 

1,000 

3 

1 

10,000 

4 

10 

15,849 

4.18 

15 

31,623 

4,5 

31 

100,000 

5 

100 

1,000,000 

6 

1,000 

10,000,000 

7 

10,000 
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Water Content (grams / lOOg Soil) 


Soil moisture suction curves for three representative mineral soils. The solid lines 
show the relationship obtained by slowly drying completely saturated soils. The dotted 
line for the clay soil is the relationship expected when a dry soil is wetted. The 
difference between the two clay Tines is due to hysteresis (modified from Brady, 1974), 


in the differences in net surface area per unit volume and the behavior 
of water at the soil solid surfaces, It was previously stressed that 
the specific surface area of soil increases as particle size decreases, 
as particle shape becomes less spherical and as clay content increases 
(thus increasing internal interlayer surface area). The remainder of 
this section will deal with the surface interaction between the soil 
solution and the soil solids. 

As was mentioned earlier, the internal structure of the clay crystal 
lattice has a charge imbalance due to the isomorphous substitution of 
lower valence cations for silicon and aluminum ions. As a result, the 
clay mineral particles have a constant negative surface charge distri- 
buted on the layer surfaces whose magnitude is dependent upon the mineralogy 
of the lattice structure. The electroneutrality of the clay particle is 
maintained by the adsorption of cations (Na , Ca , K , M g ) on the 
crystal unit layer surfaces. These exchangeable surface cations can be- 
come completely disassociated from the particle when hydrated. 

o 

In contrast, the surface charge of oxides (such as SiO or quartz) 
is determined by the incomplete terminal molecular layer at the crystal 
edge. The loss of coordinating cations at the frayed particle edge leads 
to the partial hydroxylization of the surface which, when hydrated, pro- 
duces a surface charge highly dependent upon pH. This means that for 
soil solutions very high in H-ion concentration, the surface may adsorb 
more protons than needed for neutralization and result in a positive 
surface charge (Bolt and Bruggenwert, 1976). This situation also exists 
at the broken edges of silicon tetrahedral layers in the clay minerals such 
that there exists a pH dependent potential surface charge at the frayed 
edges and a constant negative surface charge on the layer faces. 

Within the liquid layer adjacent to the charged particle surfaces, 

"the dissociated exchangeable cations tend to remain in the neighborhood 
of the spatially fixed lattice surface charge imbalances. The surface 
charge density Q of a given soil is the ratio of the cation exchange 
capacity to the specific surface, 

Q - CEC/S, meq/m 2 . (2.7) 
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The density of the surface charges determines the magnitude of an elec- 
trostatic field extending into the adjacent liquid layer. The electric 
field attracts compensating dissociated counter-ions and repels ions of 
like charge thus producing what is termed the electric double layer (Gouy, 
1910; Stern, 1924). The accumulation of counter-ions within the hydration 
shell adjacent to the solid surface is not complete, however, since the 
adsorptive force of the electric field is opposed by the diffusion of 
counterions from the surface toward the lower ionic concentration in the 
bulk solution. The opposing tendencies of entropy and electrostatic 
attraction produce an exponential decrease in counterion concentration 
as a function of distance from the soil solid surface (Figure 2.12). 

The resultant counterion atmosphere is known as the Gouy layer which 
represents a condition of minimum free energy in the system. The Gouy 
layer also consists of a distribution of co-ions (charge similar to 
that of the surface charge) whose concentration increases exponentially 
with distance from the solid surface due to electrostatic repulsion 
(van Olphen, 1963). By use of the Poisson-Boltzmann equation from 
electrostatic and diffusion theory, the distribution of positive and 
negative ion concentrations can be calculated as a function of distance 
from a solid surface of known charge density. 

(^/Cq = e (“ AE/kT ) (2.0) 

where: C-j/Cq = ratio of counterion concentration near the 

surface to that of the bulk solution 
aE = the potential energy difference of ions in the 
attractive field 

kT = the kinetic energy of the ion where k is Boltzmann's 
constant and T is the absolute temperature. 

The ion concentration of the Gouy layer is shown in Figure 2.13 
for a negatively charged particle surface. The effects of a change in 
the solute concentration of the bulk solution are shown for a frayed 
particle edge which is typified by a constant surface potential and for 
a clay unit layer face which is typified by a constant surface charge in 
Figures 2.13a and b, respectively. The lines AD and CD represent the 
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Concentration of counterions in the double layer 

Figure 2.12, Distribution and concentration of counterions in the 

electric double layer for a) condition of mini mum' Energy, 
b) condition of maximum Entropy, and c) condition of 
minimum Free Energy, actual distribution (modified from 
Bolt and Brugqenwert, 1976). 
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Figure 2.13. Charge distribution in the diffuse double layer of a negative particle at two electrolyte 
concentrations, a) constant surface potential at crystal edges and b) constant surface 
charge at clay layer faces (modified from van Olphen, 1963). 






respective counterion and co-ion concentrations in the soil solution as a 
function of distance from the particle surface. The extent of the diffuse 
double layer into the bulk soil solution is given as the distance from the 
solid surface at which the co-ion and counterion concentrations are equi- 
valent to those of the bulk solution; this is shown graphically as the rela- 
tive length of line BD. The net surface charge of the particle 
surface is equal to the area defined by CAD. In Figure 2.13a 
the Gouy layer of the frayed crystal edges is shown to be pH dependent; 
an increase in the bulk electrolyte concentration from level BD to B ' D ' 
results in the compression of the double layer towards the particle sur- 
face from distance p to y 1 and results in a net increase of the surface 
charge. In contrast, an increase of the bulk electrolyte concentration 
adjacent to a clay layer face (Figure 2.13b) causes no change in the net 
surface charge although the Gouy layer is also shown to be compressed towards 
the particle surface. It should be remembered that the effects shown in 
Figure 2.13b dominate the net behavior of the soil when clay is present 
because of the large specific surface of clay crystal layer faces having 
a constant surface charge relative to the area of crystal edges or the 
external surfaces of coarser silt and sand sized particles which have a 
constant surface potential. 

The theoretical thickness of the diffuse double layer, as defined 
above, is shown in Table 2.8 for various electrolyte concentrations of 
monovalent and divalent counterions. For normal soils, the total elec- 
trolyte concentration varies from about .01 normal at field capacity to 
roughly .1 normal at wilting point with approximately equal amounts of 
monovalent and divalent cations in solution, while saline soils are 
about .1 normal electrolyte concentration at field capacity (Bolt and 
Bruggenwert, 1976). If the diameter of a water molecule is taken as 
approximately 2.76A 0 , the thicknesses shown in Table 2.8 correspond to 
from 2 to 362 monomolecular layers of water. 

The theoretical thicknesses of the double layer given in Table 2.8 
assume that a particle surface is adjacent to an infinite volume of bulk 
solution; this condition is not usually true within the soil matrix. 

Within the soil, interparticle distances are often less than twice the 
electric double layer; Figure 2.14 compares the counterion concentration 
for such a situation with that of unimpeded Gouy layers. 
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TABLE 2.8 

Approximate "Thickness" of the Electric Double Layer 
As a Function of Electrolyte Concentration at a 
Constant Surface Potential 
(from van Olphen, 1963). 


Concentration of 
hint of op|<u.itc 
charge lo that of 
the lurlitlc, 
itiolc/litcf 



"Thickncti" of the Uotiblc layer, A 
Monovalent loot OivalcotTonT 
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Figure 2.14. The distribution of ions in the Gouy layer between two 

parallel plane micelles. Broken lines indicate unimpeded 
Gouy layers (from Childs, 1969). 



As a soil is dried by evaporation, water is removed in the vapor 
phase and ions are forced to remain in the liquid phase. Thus, as drying 
proceeds due to evapotranspirative demands, the electrolyte concentration 
of soil water tends to increase. Below a certain moisture, the voids 
between adjacent soil solid surfaces will no longer be completely filled 
with soil solution and water will consist of film layers around the soil 
particle surfaces (both external and internal). Continuing dehydration 
will cause the thickness of the water film layer to become smaller than 
the double layer thickness. Truncation of the diffuse electric double 
layer causes continued adjustment of the ionic concentrations in solution 
(Figure 2.15). A truncated double layer develops a powerful swelling 
pressure when brought into contact with additional water. This swelling 
pressure can be of sufficient magnitude in expanding lattice clays, such 
as montmoril Ionite, to close inter-aggregate pores and reduce the 
hydraulic conductivity of the soil. 

The spacing between the centers of adjacent crystal unit layers of 

dry 1:1 and 2:1 alumina silicates is approximately 7.2A 0 and 9.6A 0 , 

respectively (van Olpehn, 1963). In the case of montmoril Ionite, the 

initial contact with water produces an increased spacing on the order of 

12.5-20A 0 , depending upon the specific clay and counterions present. This 

initial interlayer swelling, at most, doubles the volume of the dry soil 

by adding 1 to 4 monomolecular water layers in a very stable configuration; 

the initial swelling is attributed to surface hydration and counterion 

hydration within the electric double layer. The forces holding water to 

the soil surface in this region are very strong; van Olphen (1963) estimates 

that the energy required to remove the last few layers of water is on the 

o 

order of 50-100 ergs/cm which corresponds to a pressure of approximately 
60,000 psi or 4,000 atmospheres. 

At crystal unit layer separations in excess of 10A° the surface hydra- 
tion energy is no longer important and subsequent swelling of clays is due 
to osmotic swelling. Osmotic swelling is caused by the double layer repul- 
sion between adjacent solid surfaces and the extent of such swelling is 
limited by van der Waals attraction between the solids and the presence 
of positively charged particle edge to negatively charged layer face 
linkages (Mitchell, 1976). Osmotic swelling reduces the truncation of 
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Figure 2.15. Concentration distribution in a truncated double. layer 
formed by removing part of the water from a system with 
an initial thickness of the water layer d-j ■ > u i while 

in equilibrium with c Q Note that j^cjix remains 

constant as no salt is’removed (from Bolt and Brugqenwert, 
1976). 




adjacent Gouy double layers to the extent permitted by the above mentioned 
inter-particle attractive forces. Osmotic swelling pressures decrease 
with increasing water content of the soil from a soil suction of about 
30 atmospheres to less than 0.1 atmosphere (van Olphen, 1963). Osmotic 
swelling can be responsible for the addition of up to 20 monomolecular 
layers of interlayer water depending upon the clay mineralogy and the 
electrolyte concentration. 

For a given particle surface of known charge density and a solution 
of known electrolyte concentration, the average electric potential can be 
computed at any point with respect to the surface from the Gouy layer model. 
The potential will have a maximum value at the surface and decrease 
exponentially with distance from the surface (Figure 2.16). An increase in 
electrolyte concentration around oxides such as quartz or clay crystal edges 
will result in a more rapid decay of electric potential with distance from 
a constant surface potential (Figure 2.16a), while the effect upon the 
constant surface charged clay layer faces is to uniformly reduce the 

electrical potential at all distances (Figure 2.16b). 

The magnitude of the electric potential at any given point in the 

soil solution determines the degree of orientation of dipolar water 
molecules within the electrostatic field maintained by the soil solid 
surfaces. Since the field strength and electric potential decrease with 
distance from the particle surface, the rigidity of water molecule orientation 
also decreases exponentially with distance (Nielson, et a!., 1972). In 
addition, the freedom of water molecules to oscillate within an impressed 
field of microwave energy is inversely proportional to the force by which 
water is rigidly oriented relative to the soil solid surface. Thus, the 
dielectric constant of water at microwave frequencies can be expected to 
increase exponentially with distance from the soil solid surface as 
electrostatic attractive forces diminish. 

The compression of the electric double layer with increasing electrolyte 
concentration (Figures 2. 13 and 2.14) is important to the present discussion 
because of its effect upon the dielectric constant of soils with low moisture 
contents and of saline soils at any moisture content. In general, the 
effect of increased electrolyte concentration will be to decrease the 
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a oca 



(a) 



(b) 

Figure 2.16. Electric-potential distribution in the diffuse double 
layer at two electrolyte concentrations, (a) Constant 
surface potential, (b) Constant surface charge (from 
van Olphen, 1963). 
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electric potential of water molecules relative to the solid surface and 
thereby increase dielectric constant. As a result, the dielectric constant 
of saline soils will be higher than that of mere neutral soils (all other 
parameters being equal), At high moisture contents of the soil where 
the volume of soil water is in excess of that needed to hydrate all of 
the diffuse electrical double layer, the excess water will have the 
dielectric characteristics of bulk water at that particular electrolyte 
concentration. 

In summation, the dielectric behavior of soil is governed by the 
dielectric constant of soil solids and the dielectric constant of soil 
water which is dependent upon the surface charge density of the solids 
at a given electrolyte concentration, the specific surface and the total 
water content of the soil. Since the surface charge densities of commonly 
occurring soil minerals are very similar in spite of large differences 
•in CEC (Table 2.9), the net dielectric behavior of soils are dominated 
by the specific surface area of soils due to particle size distribution 
and internal surface area acting in conjunction with the volume of water 
at given electrolyte concentrations. Hence, it is to be expected that 
the dielectric constant of moist soils should exhibit a pronounced 
dependence upon soil texture. 

2.2 Electrical Behavior of Soil Water 

Variance in the electrical behavior of soil water has been observed 
at microwave frequencies in measurements of the dielectric constant and 
in measurements with passive and active microwave sensor systems. These 
differences have been variously attributed to the effects of salinity, 
bulk density, and "bound" water among others. 

2.2.1 Dielectric measurements of soils 

The effects of soil texture on microwave sensor response to soil 
moisture are predicated upon the complex dielectric constant of moist 
soils. In general, the real part of the dielectric constant is about 
3.5 for the dry mineral constituents of soil and about 79.0 for pure 
water, more than an order of magnitude in difference. Dielectric 
measurements of soils have been made by numerous investigators for a 
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TAGLE 2.9 

Properties of Various Clay Minerals 
(from Baver, Gardner and Gardner, 1972). 

The Specific Surface, Cation-Exchange Capacity, and Density of 
Charge ' 


Clay mineral 

* 

Specific surface 
(rn’/K) 

Cation-exchange 

capacity 

(tne/6) 

Density of charge 
(me/m* X 10»> 

Kaolinitcs 


v* 

1 

o 

* 


0.03-0, 15t 

6-7,5) 

Hikes 


100-200’ 


0,10-0,40) 

1, 0-2.0) 

Vcrmieulites 


300-500* 


1.00-1.50) 

3.0-3. 3) 

MontmoriHonhcs 


700*800* 


0.80-1.50) 

1.1- 

•1.9) 


Ns 

a»B 

E.G. 


N, 

CPB 

lllitc-1 § 

93 

96 

91 

0.26 

2.8 

2,7 

IHitc-2 5 

132 

138 

144 

0.41 

3.1 

2.9 

Kaolinite-1 $ 

17 

21 

21 

0.043 

2.5 

2.0 

Kaolinite-2§ 

36 

36 

— 

0.050 

1.4 

1.4 

Kaolinite-3 § 

40 

9 

13 

0.030 

0.75 

3.3 

Mommoritlonhc-I § 

47 

800 

— 

0.98 

21.0 

1.22 

Mommorillonite-2§ 

49 

600 


0.98 

21.0 

1.6 

Montmorillonhe-Jlj 

101 

800 


0.99 

9.8 

1.22 


* Fripiat (1964). 

)Grim (1962), (Used with permission of McGraw-Hill Book Co,) 
) Calculated from the data in second and third columns. 

5 Greenland and Quirk (1964). 

CPB — cetyl pyridinium bromide. 

E.G. — ethylene glycol. 


Heat of Wetting, Specific Surface, Cation-Exchange Capacity, and 
Density of Charge (from Greene-Kelly, 1962). 


Clay mineral 

Heat of 
wetting 
(cal/g) 

Specific 

surface 

(m*/g) 

CEC 

(mc/g) 

Density of 
charge 

(me/m 1 X 10’) 

Kaolinite-I 

1.6 

14.8 

0.035 

2.4 

2 

1.4 

12,0 

0.023 

1,9 

3 

1.4 

11.0 

0.019 

1.9 

4 

2.1 

25.0 

0.043 

" 1.7 

Hydrous mica-1 

7.6 

150 

0,25 

1.7 

2 

4,8 

no 

0,17 

1.5 

3 

7.9 

160 

0,30 

1.9 

4 

16,5 

250 

0,43 

1.7 

Montmorillonite*l 

16.5 

690 

0.92 

1.3 

2 

17.4 

640 

0.83 

1.3 

3 

22,2 

700 

1.13 

1.6 


variety of soil textures. Results pertinent to the present discussion 
are presented below. 

The complex dielectric constant consists of a real and an imaginary 

part: 

K r * K; - JKJ (2.9) 

where: j » /T 

KJ, * the relative permittivity 

KJi * the effective conductivity 

The real and imaginary parts of the complex dielectric constant for water 
can be computed from the Debye equations as modified by Lane and Saxton 
(1952) to include ionic conductivity effects of salt solutions: 

(K, - K») 

K rw = Kw + TTJf/fE) 2 (2.10) 

and 

( K — K„) g , 

^rw “ < f / f °) T+ (f/fo) 2 + 2-rrf K q i 2 -! 1 ) 

where: KJ, W = KJ, for water 

K rw = K r for water 

K°° “ relative permittivity at frequencies much higher than 
the relaxation frequency 

f Q = the relaxation frequency 

K $ - static relative permittivity (f = o) 

= ionic conductivity of the salt solution 

K Q - free-space permittivity. 

Cihlar and Ulaby (1974) summarized many of the dielectric measurements 
of soil and noted the following as regards the general behavior of 
dielectric constant as a function of soil moisture: 
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a) Moisture Effects. At low soil moistures the real part of the 
dielectric constant increases slowly with increasing moisture to a transition 
point proportional to the clay fraction of the soil (Figure 2.17). This 
region of low sensitivity to soil moisture is attributed to the effect of 
adsorbed water on the soil particle surfaces. At moistures above the 
transition point, the real part of the dielectric constant increases at a 
much faster rate with increasing moisture. The imaginary part of the 
dielectric constant is smaller than the real part at most frequencies and 
usually increases monotonically with increasing moisture content although 
some investigators (Geiger and Williams, 1972; Newton, 1977; Wang and 
Schmugge, 1978) have noted an inflection point similar to the transition 
point defined for the real part. 

b) Bulk Density Effects. Bulk density has been shown to have a 
considerable effect on KJ, for given soils and minerals at constant gravimetric 
moisture. Cihlar and Ulaby (1974) report that studies of pumice at 10 GHz 
(Peake, et al., 1966) and of soil (Edgerton, et al . , 1971) supports the 
following relationship as shown in Figure 2.18: 

k;=(i+x )2 (2.12) 

As a result, it is necessary to express soil moisture as a volume percent 
to account for the dependence of on bulk density. Figure 2.19 illustrates 
the effect of using volumetric moisture as the soil water descriptor in 
place of gravimetric moisture. It can be seen that volumetric moisture 
units adequately account for variance in KJ, values at all measured moistures 
regardless of soil compaction. 

c) Soil Textural Effects. As previously mentioned, the primary 
effect of soil texture on the measured dielectric constant appears to be the 
value of the transition moisture. The existance of a transition moisture 

or moisture range is supported by the electric double layer model which 
contends that the first few layers of absorbed water should have a limited 
mobility and thus inhibit vibration at microwave frequencies. -There 
exists some uncertainty as to the structure of the first few layers of 
water forming a hydration shell around hydrophilic particles. If this 
water has a hexagonal grouping structure similar, but not identical, to 
that of ice (Hendricks and Jefferson, 1938; Ravina and Low, 1972), it could 
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DIELECTRIC CONSTANT 



Figure 2.47 


A comparison between the calculated dielectric constants from 
the empirical model and the measured values at 1.4 GHz. 
(From Wang and Schmugge, 1978). 


! 1 


51 





SAN ANTONIO CLAY LOAM 

• 0.03 GHz 

• 1.0 GHz 
▼ 3.8 GHz 

PUERTO RICO CLAY LOAM 
o 0.03 GHz 
□ 1.0 GHz 
A 3.8 GHz 

MOISTURE CONTENT: 0% 
SOURCE: HIPP (1974) 
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Figure 2.18 Effecf of bulk density on fhe relative dielectric constant, 
(from Cihlar and Ulaby, 1974). 





DIELECTRIC CONSTANT, REAL PART 


SOIL: OPENWOOD STREET S ILT 
FREQUENCY: 1.074 GHz 
COMPACTION MOISTURE BY 
(N/cm 2 ) VOLUME WEIGHT 


5. 74 • o 

11. 83 ▼ a 

18. 47 ■ □ 
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SOIL MOISTURE (PERCENT BY WEIGHT) 


Figure 2.19 


Change in relative dielectric constant due to soil moisture 
units used (from Cihlar and Ulaby, 1974). 



have a dielectric constant about 0.1 that of free water (Grim, 1960). 

However, there exists evidence suggesting that water within the hydration 
shell is more dense than free water and more highly ordered than ice 
(Mitchell, 1976). The double layer model also predicts that conductivity 
should increase at low moisture contents due to the presence of ion 
concentrations inversely proportional to distance from particle surfaces. 
Since the imaginary part of the dielectric constant KJ1 is proportional to 
conductivity, the magnitude of K” is proportional to the specific surface 
area of soils at a given moisture. Measurements of K" for sand and loam 
at .13 to 3 GHz (Leschanskii , et a),, 1971; Lundien, 1971) demonstrate 
that of loam (Figure 2,20) is much higher than that of sand (Figure 2.21). 
These figures also demonstrate that the difference between of loam and 
sand decreases with increasing frequency. 

There have been many theoretical mixing formulas applied to describe 
the dielectric behavior of the soil -water-air system (Table 2.10), A 
comparison between the mixing formulas reported in the literature and 
actual values measured for sand, clay loam and clay (Figure 2.22) demonstrates 
that these theoretical models are poor descriptors of the observed dielectric 
behavior and fail to account for the pronounced differences between soil 
textural classes (Wang and Schmugge, 1978). More complex empirically 
derived mixing formulas have been developed which more accurately predict the 
observed dielectric constant as a function of volumetric moisture M and 
can account for the effects of soil texture (Wobschall , 1977; Wang and 
Schmugge, 1978). Wobschall (1977) treats the soil-water-air system as 
a complex multiple mixture of various portions of the soil solution and 
soil solids which appear to be qualitatively related to the double layer 
model and soil tension. In a more explicit approach relative to soil 
physics, Wang and Schmugge (1978) separate volumetric moisture into a 
portion above the transition moisture which is treated as bulk water 
and a portion less than the transition moisture which is treated as 
adsorbed water having the dielectric constant of ice. In Figure 2.17 
the transition moisture is indicated graphically as the intersection of 
the two line segments for each soil. Transition moisture was observed 
to be independent of frequency for a given soil and was found to be 
highly correlated with an estimate of the water content at 15 bars of 


54 


_ 24 

o ; 

£ 22 
I 20 

z 


a OPENWOOD SILT LOAM, 1. 1 GHz (LUND IEN, 1971) 
a OPENWOOD SILT LOAM, 1.5 GHz (LUND|EN,1971) 
■ RICHFIELD SILT LOAM, 0.3 GHz (LUND IEN, 1966) 
o LOAM, 0. 13 GHz (LESCHANSK 1 1 et al. , 1971) 

T LOAM, 0.3 GHz (LESCHANSK 1 1 et al. , 1971) 
o LOAM, 1. 0 GHz (LESCHANSK 1 1 et al. , 1971) 

(7 LOAM, 3.0 GHz (LESCHANSK 1 1 et al. , 1971) 


0 5 10 15 20 25 30 35 40 

SOIL MOISTURE (PERCENT BY WEIGHT) 

Figure 2.20 

Relative dielectric constant vaiuos of soil as a function of gravimetric 
water content; Loam; frequency 0.13 GHz - 3,0 GHz; imaginary part 
(from Cililar and Ulaby, 1974). 
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DIELECTRIC CONSTANT 


• SAND, 0.13 GHz (LESCHANSK lletal., 1971) 
o SAND, 0.3 GHz (LESCHANSK lletal., 1971) 
x SAND, 1. 0 GHz (LESCHANSK 1 1 et al. , 1971) 
a SAND, 3. 0 GHz (LESCHANSK 1 1 et al. , 1971) 
a YUMA SAND, 1.1 GHz (LUNDIEN, 1971) 

* YUMA SAND, 1.5 GHz (LUNDIEN , 1971) 
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SOIL MOISTURE (PERCENT BY WEIGHT) 


2.21 


Relative dielectric constant values of soil as a function of 
gravimetric water content: Sand; frecuency 0.13 GHz - 
3.0 GHz; imaginary part (from Cihlar and Ulaby, 1974). 



TABLE 2.10 

The Mixing Formulas Used for Comparison with the Measured 
Dielectric Constants of Soil-Water Mixtures (from Wang and Schmugge, 1978). 


Somce 

Formula 

Bemnrks 

l. Rayleigh (1892) 

c - 1 _ r 5° ~ 1 . , £ l “ 1 
£ + 2 0 £ 0 f 2 ' ' f , + 2 

f 0 + f, » 1 

2. Butt elior (1952) 

e - £ , 

3r 0 72e 

1 

3. Brown (1956) 

e = f 0 «o + f i £ i ■ 

f o + ( r 1 

*l. Bfrchak ct. al, (1974) 

= f 0 f i 

f o + f , ” 1 

5. Bruggcman (1935) 



G. Wagner (1914) 

£ " £ i , £ o- £ , 

3c, 0 £ o + 2e, 


7, Kharadly and Jackson (1953) 

* - £ » _ { £ o - £ l 
c + 2e, u £„ + 2c, 


8. Ncrpln and Cluidnovsltii (1970) 

e, - 1 

£=1 £-1 

f * fractional volume 
oT each constituent 

9. Wiener (1910) 

r _ *o £ o^ + £ i .. _ £ i + F 

f«U + f i £ o + F 

f 0 + f, *• 1. F is an 
adjustable parameter 

10 . Pearce (1955) 

f.(l-F) 

£ = e ° = i - r,K 

-t<Fil 
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DIELECTRIC CONSTANT e' 



VOLUMETRIC WATER CONTENT (cm 3 /cm 3 ) 

Figure-2 .22 . A comparison between some mixing Formulas reported in the 
literature and the measured dielectric constants of 3 soils as 

a function of water content. Curves are labeled according 
to Table 2.10 (From Wang and Schmugge^ 1 978) . 





tension: 


W T = 0.09 + 0.59 x WP (2.13) 

3 3 

where: Wy R transition moisture, cm /cm 

WP = an estimate of moisture at the wilting point at 
15 bars of tension, cm 3 /cm 3 

WP = 0.06774 - 0.00064 x % sand +0.00478 x % clay 

From observations of the behavior of the transition moisture as 
a function of particle size distribution, Newton (1977) postulates a 
direct dependence of the complex dielectric constant on the tension 
of soil water. In lieu of experimental measurements of and as 
a function of tension, Newton (1977) assumes the moisture-tension 
characteristic of the clay loam (Figure 2.23) to be equivalent to that 
measured for the compositional ly similar Avondale clay loam. Further- 
more, assuming relative permittivity to be solely dependent upon soil 
tension in the manner given by the clay loam (Figure 2.23), the relation- 
ship between relative permittivity and tension is obtained (Figure 2.24), 
which is texture independent by definition. Combining Figures 2.23 
and 2.24 yields estimated volumetric moisture-tension characteristic 
curves for sand and clay as seen in Figure 2.25. The strong similarity 
between the families of estimated and measured moisture-tension 
characteristics. Figures 2.25 and 2.26 respectively, supports the validity 
of the postulated dependence of the dielectric constant on soil water 
tension. 

2.2.2 Observed soil textural effects on microwave sensor response to 
soil moisture 

Microwave sensor response to a target is determined by the relationship 
between the frequency, incidence angle and polarization of the sensor and 
the dielectric and geometric properties of the target. Conservation 
of energy dictates that all energy incident upon a surface must-' be absorbed, 
reflected or transmitted. 

W, = W +W + W (2.14) 

1 a p t 
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Relative Permittivity 


-■*- Miller Sand 


Clay Loam 

****■— Miller Clay 
25 r 



Volumetric Moisture, M v 

Figure 2.23 Relative Permittivity of Sand, Clay Loam and Clay as a Function of Volumetric 
Soil Moisture at 1.4 GHz (Data from Newton, 1977) 
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Figure 2.24 RELATIVE PERMITTIVITY OF A SAND, CLAY LOAM AND CLAY AS A FUNCTION OF SOIL 
TENSION AT 1,4 GHZ. (Data from Newton, 1977). 
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where: 


= incident energy 

= absorbed energy 

W = reflected energy 
P 

W «* transmitted energy 

T 

Normalizing equation 2.14by W i at a given wavelength X yields: 

1 = a(x) + p(x) + r(x) (2.15) 

where: a = absorptivity 

p = reflectivity 
t = transmittivity 
X - wavelength 

For an opaque surface, transmittance = 0 and equation (2.15) can be 
simplified to: 

1 = a(x) + P (X) • (2.16) 

Kirchhoff's Law states that absorptivity equals emissivity e. 

a(x) = e(x) (2-17) 

By substituting equation 2.17 into 2.16: 

e(x) =1 - p(x) • (2.18) 

Passive microwave radiometers measure the brightness temperature of 
the surface which is a function of the emissivity and the surface temperature, 

T b - E T g (2-19) 

where: Tg = observed brightness temperature °K 

Tg = surface temperature °K 


Active microwave sensors measure the backscattering coefficient 
a 0 of the surface which is related to the power reflection coefficient 
R at a given incidence angle 0 and the surface geometry. 

a° Q = f(R., roughness) (2,20) 

The reflectivity of a surface is governed by the geometric roughness 
relative to the sensor wavelength and the dielectric properties of the 
surface. According to the Rayleigh criterion, a surface may be 
considered electromagnetically smooth if 


h $ x/(8 cose) 


( 2 . 21 ) 


where: h = height variations above a reference plane 

0 = incidence angle (relative to nadir) 


For a smooth surface, reflectivity is equal to the power reflection 
coefficient and the square of the voltage reflection coefficient 

p = R = r 2 . (2.22) 

which is a function of the complex dielectric constant. The Fresnel 
voltage reflection coefficient r is defined separately for each 
polarization of the incident electric field by: 

r v exp(.1* r > - -" 2 cos8 i + <« 2 - s1n2e i )% (2.23) 

n coso.j + (0 - sin o^) 2 


and, 


r H exp( j4> h ) = 


O O 1 

cos e^ - (n^ - sin o.j)' 2 

— ^ ^ -j- 

cose i + (n " - sin e..) 2 


(2.24) 


where: r = voltage reflection coefficient for the vertical or 

horizontal polarization, r v and respectively. 

exp(o'i>) = the phase relation between the incident and 
reflected electric field. 

n = the index of refraction between two media. 
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Since at microwave frequencies 


n 



(2.25) 


where: “ relative permeability = 1.0 for soils. Equations (2.23) 

and (2.24) can be expressed in terms of the dielectric constant K p as: 


k r cosG - (k y - sin^o)* 5 
k r coso + (k f - sin^ep 


(2.26) 


and, 


coso - (k r - sin 2 o) 3a 

— j_ 

coso + (k - sin o) 


(2.27) 


At normal incidence angle and for a smooth surface, the above equations 
reduce to a single expression independent of polarization, 


r- 


I'fr + 1 


(2.28) 


Using the data from Figure 2.23, the relationship shown in 
Figure 2.27 is obtained for R (in dB) at an incidence angle of 0°. The 
sensitivity of R as a function of volumetric moisture M y is observed to 
be inversely proportional to clay content. Since the radar backscatter 
coefficient a 0 and the normalized antenna temperature are proportional to 
R and ( 1 - R ) respectively, it is apparent that the sensitivity of active 
and passive microwave sensor response to should also be inversely related 
to percent clay content of soil. Furthermore, if it is assumed that the 
dielectric constant is a function of tension and independent of soil texture 
as shown in Figure 2.24, then it follows that R will also be texture 
independent as a function of tension (Figure 2.28). Moreover, expressed 
in dB, R appears to be linearly related to the logarithm of tension 
(Figure 2.29). 

Experimental verification of the relationship postulated by Figure 
2.29 cannot be derived directly from prior passive or active microwave soil 
moisture experiments because soil tension was not measured as a part of the 


Reflection Coefficient at 0° (dB) 



Figure 2.27 Reflection Coefficient (dB) at 1.4 GHz, 0° as a Function of Volumetric Soil Moisture for Sana, 
Clay Loam and Clay. (Data from Newton, 1977) 


Fiqure 2.28 REFLECTION COEFFICIENT OF SAND, CLAY LOAM AND CLAY 
TENSION AT 1.4 GHZ, 0° (Data from Newton, 1977). 
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FIGURE 2.29 REFLECTION COEFFICIENT (dB) OF SAND, CLAY LOAM AND CLAY AS A FUNCTION OF SOIL 
lENSION AT 1.4 GHZ, 0°. (Data from Newton, 1977). 
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"ground truth". However, support for the above does come from the use of 
a normalized moisture M defined as: 

M n = 100M y /M v {l/3 bar tension) (2.29) 

3 

where: M y = volumetric moisture, g/cm 

M y (l/3 bar tension) - M y at 1/3 bar tension 

Passive microwave experiments near Phoenix, Arizona using a 21 and 1.55 
cm airborne radiometer showed that the sensitivity of brightness 
temperature to M y was inversely proportional to clay content as 
demonstrated in Figure 2.30 (Schmugge, 1976a, b). It was found that 
conversion of gravimetric moisture Mg to normalized moisture M R produced 
a relatively texture independent relationship. Percent of field 
capacity was estimated by Schmugge (1976a) from gravimetric moisture 
Mg by: 

M pc = 100 x Mg/FC (2.30) 

where: FC = estimated field capacity or Mg at 1/3 bar tension 

and, 

FC = 25.1 - 0.21 (56 sand) + 0.22 (56 clay) (2.31) 

While equation 2.31 is only one of many different empirical relation- 
ships (Table 2.11) derived to relate moisture at a given tension to soil 
texture, they all produce similar values of FC and the semantic ambiguity 
of "field capacity" as defined by an estimate of any laboratory determined 
water retention is overshadowed by the significance of the texture 
independent relationship displayed in Figure 2.31. The findings of Schmugge 
(1976a, b) offer additional support to the previously postulated dependence * 
of K r on soil tension since the data from Figures 2.27 to 2.29 also exhibit 
a linear dependence of R(dB) on M p (Figure 2.32). 

The above relationship is further supported by results derived from 
a series of truck-mounted radar investigations of soil moisture for both 
bare and vegetation covered test fields. In a comparison of non-vegetated 
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Brightness Temperature, Kelvins 


Figure 


T b = 310 - (2. 57 + 0. 13) X T„= 311 - (1. 45 + 0. 086) X 




€-1 cm gravimetric moisture 

2.30. Plots of 1.55 cm brightness temperature versus soil moisture for light soil (sandy 
loam and loam) and heavy soil (clay loam) for bare fields (from Schmugge* 1976). 


TABLE 2.(1 

Estimation of Field Capacity and Wilting Point from Soil 
Texture and Other Commonly Measured Soil Variables 


Definitions: Field Capacity = FC = moisture at 1/3 bar tension 

Wilting Point = WP = moisture at 15 bars tension 
Available Water Capacity = FC-WP = water available 

to plants 

Textural classes are based on USDA 1951 classification. 


I. Salter & Williams (National Vegetable Research Station, Britain) 

AWC (cm/cm) = 22.32 - 0.132 (% coarse sand) + 0.216 (% silt) 

Mean Error of Estimate = +10.0% 

Regression is based on 27 different soil textures from soils 
analyzed by the authors. 

II. Schmugge, et al . (NASA TN D-8321) 

WP (%) = 7.2 - 0.07 (% sand) + 0.24 {% clay) Multiple R = 0.945 
FC (%) =25.1 - 0.21 {% sand) + 0.22 (% clay) Multiple R = 0.904 

Regressions are based on 111 cases of 37 independent soil textures 
from soils near Phoenix, Arizona and the Rio Grande Valley in 
Texas. The soil data is derived from local Soil Conservation 
Service reports. 

III. Dobson & Ulaby 

WP (%) = 4.80 + 2.47 (% organic carbon) + 0.24 (% clay) Multiple R = 0.8 
Standard Error of Estimate = 3.66 

FC {%) = 35.29 + 3.74 {% organic carbon) - 0.30 (% sand) - 0.15 {% silt) 
Multiple R = 0.80 
Standard Error of Estimate =6.31 

Regressions are based on the surface horizon of 190 independent 
soil types from the Soil Conservation Service Laboratory Reports 

IV, Brasher, Flach and Nettleton 

FC (%) = 0.467 [% clay) + 0.356 (% fine silt) - 17.36 (bulk density 
at oven drying) + 0.039 (soil depth in cm) 

Multiple R = .842 

Standard Error of Estimate =5.7 
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Percent of 1/3 Bar Moisture, M n 

Reflection Coefficient (dB) at 1 .4 GHz, 0° as a Function of Normalized 
Volumetric Soil Moisture, Percent of 1/3 Bar Water Content, for Sand, 
Clay Loam and Clay. (Data from Newton, 1977) 


soil results from 44 observations of Miller clay and 80 observations 
of Eudora silt loam, it is evident that the sensitivity of o° at C-band 
to M y is dependent upon soil texture (Figure 2,33a) while conversion of 
M y to M pc by use of Equation (2.20) produces the texture independent 
relationship (Figure 2.33b) found by Batlivala and Ulaby (1977). Further- 
more, the value of M^q is substantiated for a 0 response to soil moisture 
under a variety of vegetation canopies including corn, soybeans, grain 
sorghum and wheat over an entire growing season from a flood plain test 
site encompassing soils ranging from sandy loam to silty clay (Ulaby, 
et al., 1979). Figure 2.34 is a scatter plot of a 0 at 4.25 GHz, HH 
polarization and 10° incidence angle as a function of Mp^ in the 0-5 cm 
soil layer ; correlation coefficient is found to be 0.917. A comparison 
of the linear regression fits of a 0 to individual crops for moisture 
expressed gravimetrically (Figure 2.35a) and as a percent of field 
capacity Mp- C (Figure 2.35b) shows a significant improvement in the 
combined correlation coefficient from 0.82 to 0.92 when Mp^ is employed 
as the soil moisture descriptor. Furthermore, the above result is not 
unique to a specific depth interval over which soil moisture is averaged; 
Mp C is demonstrably superior to either gravimetric or volumetric 
moisture over any depth interval considered (Figure 2.36). 


3.0 1977 BARE SOIL EXPERIMENT 

An experiment was conducted from August 4 to October 21, 1977 using 
the truck-mounted University of Kansas Microwave Active Spectrometer (MAS) 
1-8 GHz system on test fields located in the Kansas and Wakarusa River 
floodplains east of Lawrence, Kansas. Measurements acquired during the 
experiment are summarized in Table 3.1. The MAS 1-8 system is a calibrated 
FM-CW Radar which operated during the experiment at 8 frequencies between 
1.2 and 7.6 GHz at all three linear polarization configurations and at 
incidence angles of 10, 15 and 20° (Table 3.2). The system is fully 
described in Brunfeldt, et al (1979) and Ulaby, et al (1979); its 
calibration and accuracy are described in Stiles, et al., (1979), The 
selection of the incidence angle Q to be in the 10*-20° range was based 
upon the optimum radar soil moisture sensor parameters recommended by 
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(a) 


(b) 


Figure 2.33. Comparison of regression results of 1974 and 1975 bare soil experiments with moisture 
expressed (a) volumetrically and (b) as percent of field capacity (modified from 
Batlivala and Ulaby, 1977). 
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Backscatter Coefficient c° (dB) 


Percent' Field Capacity in the 0“5cm Soil Layer as a Function of 
Backscatter Coefficient at 4.25 GHz, HH, 10° for Corn, Milo, 
Soybean and Wheat Data Sets Combined, (Adopted from Ulaby, 





2.35 0° RESPONSE TO SOIL MOISTURE FOR THE 1975 

VEGETATION EXPERIMENTS 

0-5cm soil moisture is expressed a) gravimetrically b) as a percent of 
the estimated field capacity. o° at 4. 25 GHz, 10°, HH. 
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Figure 2.36 j Correlation coefficient as a function of depth for soil moisture 
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TABLE 3,1 

1977 Bare Soil Experiment Data Acquisition Summary 
August 4 - October 21 



1 


1 

1 


88 data sets 

4 soil types (Field capacity range 12ft - 35ft) 

All four test fields had approximately the same surface roughness 


Sensor: 


Measured characteristic 

frequency 

range 

angular 

range 

polarizattion 

Backscatter coefficient (o°) 

1-8 GHz 

10°, 15°, 20° 

HU, VV, HV 


Tarqet: 

Measured characteristic 

soil 

intervals 

measured 

number of 
sampl ing 
location 

frequency 

of 

measurement 

1 . 

Precipitation 


4 

continuous 

2. 

Temperature 


1 

continuous 

3. 

Relative humidity 


1 

continuous 

4. 

Soil moisture 

0-1 & 1-2 cm 
2-5, 5-9 & 
9-15 cm 

8 per field 
4 per field 

once per data set 
once per data set 

5, 

Soil conductivity 

0-1, 1-2 & 
2-5 cm 

1 

continuously during 
a radar data set 

6. 

Soil bulk density 

full profile 

8 

at least 3 times 
per field 

7. 

Soil texture 

full profile 

8 

once per field 

8. 

Organic matter 

full profile 

8 

once per field 

9, 

Soil water retention 
characteristics 

full profile 

1 

once per field 

10. 

Surface slope 



once per field 

11. 

Surface roughness 


1 

once per measure- 
ment day 


l 
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TABLR 3.2 

1977 HAS 1-8 System Specifications 


Pl atform Height? 26 meters 
Antennas : 

Dishes: 1.22 meter parabolic 
Feeds. Dual Polarized log periodic 
Polarizations: HH HV VV (transmit, receive) 
Incidence Anni e: 10, 15, 20 
Samples at each Angle: 20 
Radar: 

Type: FM-CW 

Modulation: Triangular 

Cent (GHz^ eqUenCieS: 1 ' 225, 1 * 625 * 2 ' 385 ’ 3 ‘ 291 > 

4.625, 5.625, 6.625, 7.625 

Frequency Sweep: 250 MHz 

Intermediate Frequency: 50 kHz 

IF Bandwidth: 20 kHz 

Transmitter power: +15 dbm 

Cal ibration : 

Internal: Delay line 

External: Luneberg lens 


For details see Brunfeldt, et al. (1979) 


Batlivala and Ulaby (1977) and Ulaby and Batlivala (1977). The scope, 
duration, test site, data acquisition procedures and data products of 
the 1977 Bare Soil Experiment are fully documented in Dobson (1979). 

3.1 Experiment Description 

The purpose of the 1977 Bare Soil Experiment was to verify that the 
1-8 GHz radar response to gravimetric and volumetric soil moisture is 
dependent upon soil texture and provide a comprehensive data base capable 
of producing more versatile soil moisture estimation algorithms that 
minimize the error component related to soil texture. 

Prior research has established the following list of target parameters 
as having probable significance in determining the 1-8 GHz radar back- 
scattering characteristics of the soil surface: 

a) soil moisture 

b) soil surface roughness 

1. surface gradient relative to incidence angle of sensor 

2. macroroughness as defined by the periodic ridge-furrow 
patterns induced by agricultural implements 

3. microroughness as defined by the size distribution of 
surface aggregates such as clods and soil peds 

c) soil bulk density 

d) soil composition 

1. texture as defined by the particle size distribution 
of the soil mineral fraction 

2. organic matter content 

3. soil salinity 

e) vegetation canopy 

1 . type 

2. height 

3. density of biomass 

4. moisture content 

5. plant geometry 

The success of the experiment in regard to stated purpose was 
predicated upon the selection of test fields which were internally homo- 
geneous and similar with respect to all target parameters, except soil 
texture and soil moisture, for the duration of the experiment. Test 
fields were selected and groomed such that they: 
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1 ) were devoid of a vegetation canopy, 

2) had surface gradients of less than .5%, 

3 ) were smooth with respect to mncroroughness, and 

4 ) had similar microroughness. 

These field conditions were attained by repeated disking of each field 
followed by a one to three week stabilization period prior to radar 
observation during which time events of saturating rain and strong wind 
further smoothed the soil surface and acted to stabilize the soil structure 
of the disturbed surface layer. Table 3.3 summarizes the initial field 
preparations and final microroughness (quantified as root mean square 
height) . 


TABLE 3.3 

field Preparations for 1977 Bare Soil Experiment 


Field 

Preparation 

Date 

Number of 
Times Disked 

Date of 

1 st observation 

Final RMS 
Height (cm) 

#1 

7/23/77 

4 

8/4/77 

1.3 


7/23/77 

8 

8/15/77 

0.8 

t;3 

8/20/77 

9 

9/3/77 

1.0 

H 

7/23/77 

11 

8/15/77 

0.7 


The microroughness of fields observed during the 1977 Bare Soil 
Experiment compare with only the smoothest of bare soil fields previously 
examined by the MAS 1-8 system in 1974 and 1975 (Table 3.4). 

Fields of sandy loam, silty clay loam and silty clay were examined 
over a broad range of moisture conditions. Importantly, test fields for 
each soil texture were found to be internally homogeneous with respect 
to soil texture and bulk density and will be discussed in more detail 
in the following section. 
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TABLE 3.4 

Comparison of Mean RMS Heights for Test Fields 
for 1974, 1975 and 1977 Bare Soil Experiments 


Year 

Field 

Designation 

Mean RMS Height (cm) 

1974 

R 

4.3 


M 

2.6 


S 

0.88 

1975 

1 

4.1 


2 

2.2 


3 

3.0 


4 

1.8 


5 

1.1 

1977 

1 

1.32 


2 

0.79 


3 

1.03 


4 

0.70 


The experiment was conducted during an abnormally wet fall and 
resulted in surface soil moisture conditions skewed toward conditions 
wetter than those found during either the 1974 or the 1975 experiments 
(BatliVala and Cihlar, 1975; Batlivala and Dobson, 1976). Table 3.5 
shows the mean gravimetric soil moistures from each of the three experi- 
ments and their associated standard deviations. The prevalence of wet 
conditions during the MAS 1-8 1977 experiment signifies that it will 
augment the moisture range of the combined 1974 and 1975 bare soil data 
base in the derivation of soil moisture algorithms. The larger standard 
deviation about the mean 1977 moistures reflects the variety of textures 
included in the experiment in contrast to the single textures examined 
in both the 1974 and 1975 experiments. 


mb 
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TABLE 3,5 

Mean Gravimetric Moisture Conditions for the 
1974, 1975 and 1977 Bare Soil Experiments 


YEAR 


Depth (cm) 


1974 

Mean S.D. 


1975 

Mean I S.D 


1977 

Mean I S.D 


0-1 

0-2 

0-5 

0-9 

0-15 

Number of 
Textures 

Number of 
Data Sets 


13.53 

15.74 

19.30 

20.20 

21.10 


7.63 

6.86 

5.70 

5.23 

4.65 


44 


16.18 

18.19 

19.97 

20.44 

20.41 


8.37 

5.63 

3.88 

2.78 

2.27 


80 


21.85 

24.98 

29.23 

31.02 

30.83 


13.25 

11.91 

10.67 

9.88 

9.06 
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3.2 Target Characteristics and Homogeneity 

Extensive measurements of soil texture, bulk density and surface 
roughness were made during the experiment to establish the homogeneity 
of these parameters within each test field. The data acquisition procedures 
and results are fully treated in Dobson (1979) and will only be summarized 
in the following sections. 

3.2.1 Soil texture 

A total of 28 samples were collected from each test field at the 
soil moisture sampling locations and depths indicated on Figure 3.1. Soil 
textural separates were defined by the hydrometer method and organic matter 
determinations were made using a variant of the Walkley-Black method of 
oxidation with chromic acid, (Allison, 1965). The values of organic matter 
content are estimated at 1.724 times the organic carbon content. A 
summary of the mean 0-15 cm soil texture and organic matter content from 
each field is given in Table 3.6 
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Figure 3.1 1977 Bare Soil Moisture Experiment: Field Layout and Sampling Grid 
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Soil Sampling Locations: 


0 - 1 cm and 1 - 2 cm samples only (locations 2 , 4, 5, and 8) 

0 -..1, 1 “ 2 , 2 - 5, 5 ~ 9, and 9 ~ 1 15 cm samples (locations 1, 3, 6, and 7) 

Location of soil conductivity sensors for 0 “ 1 , 1 ” 2, and 2 - 5 cm depth intervals 
Indicates maximum area illuminated at 1.1 GHz 

Dashed lines indicate the center of illumination by the MAS 1 ** 8 at various angles of incidence 
Elevation beamwidths vary with frequency from 12° at 1 . 1 GHz to 2° at 7.25 GHz. 





TABLE 3.6 

Summary of Organic Matter Content and Soil Texture of Test Fields 


Field 

0-15 cm Mean 

0-15 cm Mean Soil 

Texture 

0 

urycUiiL- i la LLuf uun lcii u 

(% by weight of dry soil) 

ESI 

• % Silt 

% Clay 

Soil Textural Class 

#1 

0.7 

55.0 

42.2 

2.8 

Sandy loam 

#2 

3.1 

9.1 

53.9 

37.0 

Silty clay loam 

#3 

3.9 

11.5 

52.2 

36.2 

Silty clay loam 

#4 

3.8 

' 

9.0 

45.1 

45.8 

Silty clay 















Within each field, the textural separates within each of the sampled 
layers (0-1, 1-2, 2-5, 5-9 and 9-15 cm) were treated as a sample population. 
The sample means and standard deviations are presented in Table 3.7. An 
indication of textural homogeneity within a field is given by the low 
standard deviation about the sample means which is typically less than 3 % 
with the notable exception of the top one cm layer of field 2 where the 
standard deviations exceed 6?!. 

Skewness shows the degree to which a distribution approximates a 
normal curve since it measures deviations from symmetry. The skewness 
of a normal distribution is zero, while a positive or negative skewness 
indicates that cases are clustered below or above the mean respectively. 
Kurtosis measures the relative peakedness of the curve defined by the 
distribution of cases. The kurtosis of a normal distribution is zero, 
while a positive or negative kurtosis indicates that the distribution 
is more peaked or flatter than a normal distribution. Calculation of 
the skewness and kurtosis of all textural sample distributions indicates 
approximately normal but ^latter distributions. In only one case, 
field 1, 1-2 cm clay, is kurtosis found to exceed 0. Values fall within 
the following ranges: 

* 

-0.737 s skewness $ 1.268 and 
-2.428 $ kurtosis $ .301 

The above findings are consistent with an assumption of field homogeneity 
with respect to soil texture with the exception of the surface layer of 
field 2 where -considerable variance was observed over short distances 
especially for sand and clay fractions (up to 14% change between adjacent 
sample locations in Figure 3.1). 

Fields 2 and 3 have very similar mean textures; both are silty 
clay loams by the USDA classification. This conclusion is supported 
by t-tests run between the field sample distributions for each ..depth 
interval. The validity of the t-tests is based upon assumptions that 
the sample distributions are both normally distributed and have equal 
variances. The first assumption is validated by the skewness and kurtosis 
measurements of each sample distribution, and the validi J ,y of the second 
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Mean Soil 


a. Field 1 



DEPTH (CM) 

" SAND 


MEAN 

S.D. 

0-1 

55.00 

2.56 

1-2 

55,33 

3.31 

2-5 

56.00 

1.73 

5-9 

53.50 

2.58 

9-15 

55.50 

1.15 

0-15 

Weighted Mean 

55.03 


b, Field 2 
DEPTH (CM) 

% SAND 



MEAN 

S.D. 


1 

0-1 

11.26 

6,14 

1-2 

10.48 

2.45 


2-5 

9.63 

2.71 

j 

5-9 

9.78 

0.51 


9-15 

7.71 

3,25 

\ 

0-15 

9.08 



* Weighted Hear. 


c. Field 3 

DEPTH (CM) X SAND 



MEAN 

S.D. 

0-1 

10.87 

1.93 

1-2 

11.92 

1.31 

2-5 

12,95 

1.00 

5-9 

11.15 

1.92 

9-15 

11.11 

1,93 

0-15 

11.53 



Weighted Mean 


d. Field A 

DEPTH (CM) X SAND 



MEAN 

S.D. 

0-1 

8,87 

2.71 

1-2 

8,31 

2.11 

2-5 

10,42 

2.05 

5-9 

9.70 

2,57 

9-15 

8.06 

2.22 

0-15 

9.04 



Weighted Mean 



1.7 

at Each Depth 


% SILT * CLAY 


MEAN 

S.D. 

MEAN 

S.D. 

43.00 

2.20 

2.00 

1.20 

41.38 

3.23 

2.75 

1.67 

41.75 

1.89 

2.25 

0.50 

43.00 

3.37 

3.50 

2.16 

41.75 

1.26 

2.75 

0.96 

42,18 


2.80 



% SILT % CLAY 


MEAN 

S.D. 

MEAN 

S.D, 

55.53 

6.84 

33.25 

6.14 

56.28 

3.23 

33.25 

2.78 

54.43 

4.38 

35.90 

1.73 

55.25 

1,26 

34.98 

1.67 

52.11 

5.16 

40.21 

4.35 

53.92 


37,02 



X SILT * CLAY 


MEAN 

53.33 

51.16 

52.53 

52.70 

51,79 

52.24 

S.D. 

2.03 

0.65 

3.88 

3.85 

2.72 

MEAN 

35.75 

36.90 

34,50 

36.15 

37.11 

36.23 

S.D, 

1.68 

1.75 

3.34 

2.35 
3.80 

% SILT 

X CLAY 

MEAN 

S.D. 

MEAN 

S.D. 

46.90 

4.67 

44,21 

4.50 

46.31 

2.33 

45.36 

1,99 

46.31 

5.66 

43.27 

5,14 

42.35 

2,79 

47,93 

5,14 

45,88 

3,74 

46.05 

4.50 

45.12 


45.83 
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assumption v/as checked using variance ratio or F-tests. F- tests rejected 
the hypothesis of equal variance in 40% of the between field comparisons 
at a 0.10 significance level where in excess of 60% of the rejections 
involved field 2 which tended to have a higher variance. The results 
of the t- tests are summarized in Table 3.8 for the top cm and the 0-15 
cm layer; other depth intervals do not deviate significantly from these 
results, 

In summary, all fields except field 2 were found to be internally 
homogeneous with respect to soil texture and fields 2 and 3 were found to 
represent the same textural population. Cumulative plots of the particle 
size distribution from a representative sample of each field are shown in 
Figure 3.2 and are summarized as a histogram in Figure 3.3. In addition, 
all fields were found to have comparable organic matter contents of less 
than 4.0% by weight. 

3.2,2 Surface Roughness 

A total of 116 roughness profiles were obtained from the four test 
fields using a five foot aluminum panel. A subset of 24 profiles were 
digitized along the transects at a maximal spacing of 0.15 cm. The 
profiles were chosen to represent the initial, median and final roughness 
conditions of each field both parallel and perpendicular to the radar 
look direction. Plots of the initial roughness conditions are shown in 
Figure 3.4. The root mean square variations in surface height of all 
digitized profiles are given in Table 3.9. These measurements confirm 
the field observation that all fields became increasingly smooth as a 
function of time, especially fields 2 and 3. The comparatively large 
RMS roughness of field 1 is primarily the result of residual small scale 
macroroughness produced by disking while the persistent comparatively 
small RMS roughness of field 4 was due to the presence of surface crumbs 
and small aggregates which proved to be very cohesive and resistant to 

crumbling. The marked change in the roughness of fields 2 and 3 as a 
function of time is largely due to their intermediate textural composition 

(compared to fields 1 and 4); initial clods and aggregates possessed only 

weak structure and thus disintegrated with time. 

Analysis of the radar data assumes that the mean surface roughness 
of all fields are equivalent. This assumption was tested with t-tests and 
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TABLE 3.8 

Results of t-tests of Soil Textural Distributions 
Between Fields. Two-tailed Test, Level of Significance = 0.01 


Depth Textural Fields which cannot be rejected 

Interval Separate as having unequal mean textures. 


0-1 cm 

sand 


silt 


clay 

0-15 cm 

sand 


silt 


clay 






Cumulative Percent 


Figure 





02 0.05 0.1 0.2 0.5 1.0 

ij Coarse ( Very j Fine 'Medium} Coarse 
! Silt i Fine I Sand 1 Sand ! Sand 
Sand 

5 Diameter (mm) 

2 distribution for soils examined during the 1977 




Figure 3.4 

Initial Soil Roughness Conditions on Bare Soil 1977 Fields 1 to 4. 
Plots of Detrended Soil Surface vs. Surface Height. Transects are 
Parallel to Radar Look Direction. 


a) Field 1 Date: 8/15/77 RMS Height: 1. 66cm 



Transect Length (cm) 


b) Field 2 Date: 8/15/77 RMS Height: 1.49cm 
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c) Field 3 Date: 9/3/77 RMS Height: 1,75cm 
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m.mrw, 


1977 Bare Siil 


Field Number 

1 


Direction of Transect 
(Relative to radar 
look direction) 

• 

Parallel 

Perpendicular 

Para! 1 el 

DATE; 8/15 

1.66 

1.33 

1.49 

8/22 



0.42 

9/2 

1.15 

0.85 


9/3 




9/6 




9/20 

1.17 

1.04 


9/24 




9/28 



0.47 

9/30 




MEAN 

1.32 

1.07 

0.79 


i 

i 


:LE 3.9 


RilS Surface Roughness (cm) 

2 

Perpendicular] Parallel] Perpendicular] Parallel] Perpendicular 

0 . 9 * 

0.61 

Q.5E 
0.7C 







variance ratio tests. The variance ratio test results indicated that the 
variances about the four field means could not be considered unequal at 
the 0.1 level of significance, Two-tailed t-tests of the sample means at 
the 0.01 significance level indicated a high probability that only fields 
1 and 4 are different with respect to RMS height* all other field pairs 
were found to be similar at the stated significance level. 

Inspection of the mean values shown in Table 3.9 indicates that, on 
the basis of RMS surface roughness, the four test fields cluster into two 
groups: 

a) fields 1 and 3 with a net mean RMS height in all look directions 
of 1.19 cm, and 

b) fields 2 and 4 with a net mean RMS height in all look directions 
of 0.71 cm. 

3.2.3 Soil Bulk Density 

Soil bulk density is not easily measured for the thin soil layers 
sampled during the bare soil experiment. Because of the time involved 
in the careful extraction and delicate handling of an accurate bulk density 
sample, adequate bulk density sampling was adjudged to consist of a minimum 
of three samples at each grid location and depth shown on Figure 3.1. Bulk 
density samples were taken when field moisture profiles permitted minimal 
soil compaction during sampling. The drainage characteristics of field 1 
were such that bulk density profiles could be taken in conjunction with 
soil moisture samples during 60% of the data sets, while the other three 
fields could only be sampled for bulk density three to six times during 
the course of the experiment. Table 3.10 shows the number of bulk density 
samples acquired at each location and depth during the course of the experi- 
ment. The high percentage of clay present in Field 4 and the saturated 
soil profiles common to this field made accurate bulk density sampling 
of the surface layers nearly impossible during the experiment and is 
responsible for the low sample size presented in Table 3.10, 

Two methods were employed to collect bulk density samples.. The 
standard method was to extract core samples of known volume and then dry 
the soil to compute bulk density. When the above method proved to be 
unacceptable due to soil compaction by the coring tool, bulk density was 


TABLE 3.10 


* 


Number of Bulk Density Samples 
Collected at Each Sample Grid Location 


Field 

Depth 

(cm) 

J 

2 

3 

Locations 
4 5 

6 

7 

8 

Total 

1 

0-1 

14 

14 

14 

13 

12 

15 

14 

13 
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1-2 

15 

14 

14 

13 

13 

15 

14 

13 

111 


2-5 

15 


14 



14 

14 


57 


5-9 

12 


10 



9 

11 


42 


9-15 

11 


10 



9 

11 


41 

2 

0-1 

5 

5 

5 

5 

4 

5 

4 

3 

36 


1-2 

5 

5 

6 

6 

5 

6 

3 

2 

38 


2-5 

9 


8 



8 

7 


32 


5-9 

6 


6 



5 

4 


21 


9-15 

7 


6 


* 

5 

5 


23 

3 

0-1 

3 

4 

5 

4 

3 

3 

3 

3 

28 


1-2 

3 

4 

5 

4 

3 

3 

3 

3 

28 


2-5 

5 


5 



6 

6 


22 


5-9 

4 


6 



5 

5 


20 


9-15 

4 


6 



5 

5 


20 

4 

0-1 

1 

1 

2 

1 

1 

1 

2 

1 

10 


1-2 

1 

2 

1 

1 

1 

1 

2 

1 

10 


2-5 

3 


1 



— 

3 


7 


5-9 

7 


4 



5 

8 


24 


9-15 

T 


5 



5 

8 


25 
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determined from the liquid volume displacement of saran coated soil clods 
extracted from the locations and depths given in Figure 3.1. 

A summary of the depth weighted mean bulk densities of each field is 
presented in Table 3,11a. The standard deviation associated with each of 
these field means is given in Table 3.11b. These statistics are derived 
from the distribution of the mean bulk densities at each location within 
a field. 

Initial analysis of the above bulk density data indicates that the 
four test fields are relatively homogeneous. These results validate the 
assumption that the soil fabric of the individual test fields had stabilized 
prior to the beginning of the experiment. The actual density values of 
each field are highly deperdent upon soil texture, and hence, a low intra- 
field variance in bulk density supports the contention of low intrafield 
variance in soil texture. 

However, a more complex picture emerges when consideration is given 
to the distribution of all samples given in Table 3.10 for a particular 
field and depth. Table 3.12 presents the sample distributions of bulk 
density in the 1-2 cm layer of fields 1 and 2. The tendency toward a 
bimodal distribution of density as seen in Table 3.12b for field 2 is 
typical of fields 2, 3 and 4. The characteristics of all such sample 
distributions are summarized by statistics for the mean density, 90% 
confidence interval, skewness and kurtosis given in Table 3.13. The 

■a 

large 90% confidence interval around the mean, commonly .2 to .3 g/cm , 
is due to the combined effects of sampling error, locational variance, 
density changes in the surface layer associated with the development of 
soil structure through time and the functional dependence of density 
on moisture in fields with a high clay content such as fields 2, 3 or 
4. The values of skewness and kurtosis indicate that densities are 
normally distributed about the means and the distributions are generally 
less peaked than a normal distribution; this is due primarily to the 
fact that, with the exception of field 1, density was sampled oyer a very 
small range of gravimetric moisture especially in the surface layers. The 
negative values of skewness obtained primarily for the deeper depth interval 
of the high clay fields is indicative of the correlation between density 
and moisture which was found to be as high as -0.7 for these layers. 
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TABLE 3.11 

Comparison of Mean Bulk Densities from 
Fields 1 to 4 of the 1977 Bare Soil Experiment 


A. Weighted Mean Bulk Density 0-1, 0-2, 0-5, 0-9 and 0-13 cm 
as computed from measurements at each location. 


Depth Interval 
(cm) 


Field Designation 


0-1 

0-2 

0-5 

0-9 

0-15 


43 
42 
1.27 
1.32 
1.43 


20 
24 
1.11 
1.21 
1.30 


30 
34 
1.17 
1.27 
1.33 


B. Standard Deviation of Location Means around the field 
bulk density at each depth. 


Depth Interval 
(cm) 


Field Designation 


0.04 

0.05 

0.04 

o.o; 

0.03 

0.06 

0.05 

0.1 

0.01 

0.03 

0.04 

0.0 

0.03 

0.02 

0.02 

0.0 

0.02 

0.02 

0.01 

0.0 
















Mean B 1. 393 


Variance = 0.021 


Median • 1. 390 St. Dev. » 0. 145 


Mode ■ 1. 376 


Skewness = 0. 076 



(a) Fieid 1, Sandy Loam 1-2 cm Layer 


10 


f 

i 



Mean - 1. 207 Variance - 0. 030 
Median = 1.203 St. Dev. =0.172 
Mode = Not Unique Skewness = 0. 088 
Maximum = 1. 533 Kurtosis = - 1. 289 
Minimum = 0.932 No. of Values = 38 

Range = 0.601 



0.8 1.0 1.2 1.4 1.6 

Bulk Density (g/cm 3 ) 

(b) Fieid 2, Silty Clay Loam 1-2 cm Layer 


TABLE 3.12 

Sampling Distribution of Bulk Density in the 1-2 cm Layer for a) Field 1, 
sandy loan) and b) Field 2, silty clay loam. 
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TABLE 3.13 

Skewness, Kurtosis and 90S Confidence Limits 
Around Mean Bulk Density at Each Depth Interval 


Field 

Depth 

Interval (cm) 

Mean - 
Density (g/cm ) 

903 Confidence 3 
Interval {+/-g/cm ) 

Skewness 

Kurtosis 

1 

0-1 

1.435 

.29 

.636 

.574 


1-2 

1.393 

.24 

.076 

-.370 


2-5 

1.137 

. ,25 

.446 

.640 


5-9 

1.257 

.15 

.042 

-.626 


9-15 

1.433 

.17 

.350 

-.231 

2 

0-1 

1,219 

.26 

.172 

-.260 


1-2 

1.207 

.28 

.087 

-1.289 


2-5 

.846 

.18 

.155 

-.636 


5-9 

1.267 

.16, 

.376 

-.421 


9-15 

1.462 

.11 

-.049 

-.905 

3 

0-1 

1.209 

.30 

-.068 

-.878 


1-2 

1.260 

.24 

.877 

.478 


2-5 

1.008 

.21 

.087 

-.778 


5-9 

1,343 

.20 

-.318 

-1.343 


9-15 

1.439 

.14 

.157 

-.723 

4 

0-1 

1.291 

.19 

.154 

-.945 


1-2 

1.348 

.27 

.220 

-1.806 


2-5 

1.053 

.23 

.138 

-1.403 


5-9 

1.370 

.18 

-.958 

.446 


9-15 

1.423 

.14 

-.570 

.944 
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Since most of the bulk density samples for fields 2, 3 and 4 were 
obtained over a relatively small range of gravimetric soil moisture, 
possible shrink-swell effects upon bulk density as a function of moisture 
are not accounted for by use of a mean field bulk density. Thus, the use 
of a mea,' bulk density as an assumed constant for each field and depth 
Interval may introduce a bias into the volumetric moisture values. A 
simple physical model to account for shrink-swell effects on bulk density 
is presented in Appendix A, 

3.2.4 Soil Moisture 

Soil samples were extracted from the sampling grid locations shown in 
Figure 3.1 for the determination of gravimetric moisture in conjunction 
with each radar data set. The soils were dried to equilibrium using a 
microwave oven as documented in Dobson, 1979. 

For each data set, a mean moisture and standard deviation about the 
mean were calculated for each depth interval (0-1, 1-2, 2-5, 5-9 and 9-15 cm) 
which could then be ratioed as the coefficient of variation. These three 
statistics are summarized for the 0-T cm layer of each field in Figure 3.5. 

The wide range of surface moisture observed during the experiment was due 
to repeated heavy rainfall incidents (3 in excess of 2 inches) which allowed 
surface soil moisture to vary from the wilting point at 15 bars of tension 
to saturation (in fields 2, 3 and 4). The large standard deviations around 
the mean moisture found predominantly in fields 2, 3 and 4 tended to occur 
shortly after precipitation events during the initial drying stages of 
the surface. Surface variance in soil moisture was heightened by the 
intersection of the water table with the soil surface in fields 2 and 4 
which produced partial flooding within topogrhphically low portions of the 
test fields. A comparison of the 0-1 cm moisture distributions with those 
at the 2-5 and 9-15 cm depths (Figures 3.5, 3.6 and 3.7, respectively) indi- 
cates that the large temporal variance in surface moisture approaches a field 
constant as depth increases. Also, the standard deviation of sampled soil 
moisture is observed to decrease with depth. An examination of,these 
figures clearly indicates that soil moisture falls into three sample 
populations which become more distinct with increasing depth where field 1 
(sandy loam) is driest, fields 2 and 3 (silty clay loams) are wet and 
field 4 (silty clay) is slightly but significantly wetter. This is to be 

103 ' 


FREQUENCY FREQUENCY FREQUENCY 

(# OF DATA SETS) (# OF DATA SETS) <# OF DATA SETS) 


HISTOGRAMS OF 0-1 cm MEAN SOIL MOISTURE, STANDARD 
DEVIATION AND STANDARD DEVIATION TQ MEAN RATIO 

. ■FIELD 1 

r a FIELD 2 

a FIELD 3 

V b FIELD 4 




MEAN SOIL MOISTURE (% DRY WEIGHT) 



2 4 6 8 10 12 14 16 18 20 

STANDARD DEVIATION OF MEAN SOIL MOISTURE 





FREQUENCY FREQUENCY FREQUENCY 

OF DATA SETS) {# OF DATA SETS) {# OF DATA SETS) 


Figure 3.6 HISTOGRAMS OF 2-5cm MEAN SOIL MOISTURE, STANDARD 
DEVIATION AND STANDARD DEVIATION TO MEAN RATIO 




STANDARD DEVIATION OF MEAN SOIL MOISTURE 
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expected and is a function of the drainage characteristics and textural 
composition of each field. Field 1 drained readily due to its height above 
a roadside drainage ditch and the high hydraulic conductivity of the sandy 
soil while fields 2, 3 and 4 had poor drainage and high clay contents 
which resulted in greater water retention. 

The soil moisture histories of each field are given in Figures 3.8 
to 3,11 to show graphically the change in the moisture profile of each 
field as a function of time and are annotated as to the presence of flooding. 

3, 3 AnguJI or Response of Radar Backseatter 

The radar buckscattering coefficient o° was measured at three angles 
of incidence. er° is plotted as a function of incidence angle o for a range 
of soil moistures at 1.6, 4.6 and 7.6 GHz for fields 1 to 4 in Figures 3.12 
to 3.15 respectively. In general , the rate of decrease in o° as a function 

of 0 indicates that all of the fields are relatively smooth within the 1-8 
GHz frequency range and that all fields have similar roughness. The Ray- 
leigh criterion, equation (2.21), establishes the values for a smooth sur- 
face given in Table 3.14 at 1.6, 4.6 and 7.6 GHz at 10° and 20° incidence 
angles. For fields with RMS heights less than these values the soil sur- 
face should appear electromagnet!’ cal ly smooth and exhibit a large drop in 
ct° with increasing incidence angle o, while the o° vs. o curve should be- 
come increasingly flat for RMS heights greater than the criterion values in 
Table 3.14, The sensitivity of o° to o for a given field and moisture 
condition tends to decrease as frequency increases especially for fields 
1 and 3 which were slightly rougher than fields 2 and 4 in terms of RMS 
height', this is to be expected from the Rayleigh criterion for surface 
roughness , 

Figures 3.12 to 3.15 indicate a clear dependence of a 0 on soil 
moisture within each field at all angles and frequencies. However, as 
soil moisture conditions approach and exceed saturation, there is a 
marked drop of a° to values comparable to those of moderately wet conditions 
(Figures 3,14 and 3.15) especially at 7.6 GHz. Saturation moisture Mg as 
defined by soil porosity is approximately 32, 44, 45 and 39% for the 0-1 
cm layer of fields 1 to 4, respectively. This phenomenon can be explained 
as a roughness effect; at saturation the soil appears as a smooth surface 
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DATA SET (IN CHRONOLOGIC ORDER) 

Figure 3.10 TEMPORAL SOIL MOISTURE FLUX, FIELD 3: Silty clay 
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Figure 3.12. Angular response of o° at 1.6, 4.6 and 7.6 GHz to, a range 
of surface moisture conditions in Field 1, sandy loam. 
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Figure 3.13, Angular response of 0 0 at 1.6, 4.6 and 7.6 GHz to' a range 
of surface moisture conditions in Field 2, silty clay loam. 
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Figure 3.14. Angular response of o° at 1.6, 4.6 and 7.6 GHz to a range 
of surface moisture conditions in Field 3, silty clay loam. 
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Figure 3.15. Angular response of o° at 1.6, 4.6 and 7.6 GHz to t a range 
of surface moisture conditions in Field 4, silty clay. 
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TABLE 3.14 

The Rayleigh Criterion for Electromagnetic 
ally Smooth Surfaces, A smooth surface 
has variations relative to a reference 
plane < to the values shown for 1.625, 

L 4.625 and 7.625 GHz at 10° and 20°. 

4i 

!'■ 


Criterion Values (cm) 


it 

! ; 


20 


0 


i 


Frequency (GHz) 


10 ° 


and reflection is more specular, hence the backscattered portion of the 
signal falls more dramatically as e Increases. At nadir it is expected 
that a 0 of saturated soil would be higher than o° of unsaturated soil; 
unfortunately 0° data was not acquired during the course of the experiment. 
Other polarization configurations display angular response characteristics 
similar to those shown for HH polarization, 

3.4 Frequency Response of Radar Backscatter 

Radar backscatter was measured at eight center frequencies between 
1 and 8 GHz. The frequency responses of the like polarizations, HH and 
VV, were found to be almost Identical while the response of the cross- 
polarized antenna configuration was difficult to characterize because a 
significant number of the measurements were below the noise floor of the 
system. 

The HH polarized frequency response of a 0 is shown in Figures 3,16 
to 3.19 for each field over a range of 0-1 cm moisture conditions. For 
soil moistures below saturation, o° is observed to increase monotonically 
as a function of frequency if small fluctuations on the order of +/- 1 dB 
due to signal fading are ignored.' This behavior was expected based upon 
similar results from data sets acquired in 1974 and 1975 for non-vegetated 
fields of Miller clay and Eudora silt loam respectively using the MAS system, 
0 ° is observed to be proportional to soil moistures below saturation at 
all frequencies and incidence angles. However, soil moistures above 
saturation (Figures 3.18 and 3.19) produce a drop in a 0 relative to wet 
but unsaturated moisture conditions at all frequencies and angles. As 
previously stated, this is felt to be caused by more specular reflection 
from the saturated soil relative to the unsaturated state. Furthermore, 
when moisture conditions approach saturation, a 0 is observed to loose 
its monotonia dependence on frequency and tends to oscillate in excess 
of theoretical fading considerations. This non-mo.notonic behavior 
may be a resonance effect dependent upon frequency, incidence, 
angle and the soil moisture profile as a function of depth. Additional 
experimentation is necessary to adequately account for this behavior in terms 
of moisture profile shape and moisture layering phenomena as they affect 
depth penetration and multiple reflections of the incident beam within 
the 1-8 GHz frequency region. 
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Figure 3.16. Frequency response of o°, HH polarization at 10°, 15° and 
20° for a range of surface moisture conditions in Field 1, 
sandy loam. 
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Figure 3.18. Frequency response of o", HH polarization at 10°, 15° and 
20° for a range of surface moisture conditions in Field 3, 
silty clay loam. 
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Figure 3.19. Frequency response of o°, HH polarization at 10°, 15°, and 
?.0 1 ' for a range of surface moisture conditions in Field 4, 
silty clay. 


Also, field 2 tended to exhibit highly irregular frequency response 
at all intermediate moisture conditions although this is not readily 
apparent from examination of Figure 3.17. This was felt to be caused by 
variance in surface moisture, density and texture within the field and 
will be treated in Section 3.7. 

3.5 Radar Backscatter Dependence on Gravimetric Soil Moisture 

This section will consider the radar response to gravimetric soil 
moisture within each field separately; other soil moisture indicators 
such as volumetric moisture, normalized moisture and soil tension will be 
treated in the next chapter. 

In the previous two sections a strong dependence of a 0 on soil 
moisture was observed at all frequencies and angles. The linear dependence 
of o° at 4.6 GHz, 10° incidence angle and HH polarization on gravimetric 
moisture Mg is shown in Figure 3.20. The linear correlation coefficient 
is quite good for fields 1 and 4, .898 and .869, respectively, and also 
for field 3 at moistures below saturation, r = .844. Both fields 3 and 
4 exhibit a leveling or a drop in cr° as moisture approaches and exceeds 
saturation; as previously discussed, this is due to an increase in specular 
reflection at extremely high moisture conditions. 

In contrast, field 2 exhibits poor sensor response to mean soil 
moisture with much scatter apparent in the data and a low resultant 
correlation coefficient of approximately 0.4. Because of the similarity 
between fields 2 and 3 in terms of soil texture, both are silty clay 
loams, sensor response was also expected to be similar. Poor sensor 
response to field 2 when compared to the other test fields is relatively 
independent of the sensor parameters (frequency and angle) and of the 
depth interval of soil considered. Table 3.15 presents the linear 
correlation between a 0 at 6 frequencies at 10° and 20° HH polarization 
and the 0-1, 0-2 and 0-5 cm mean gravimetric soil moisture of each field. 
Sensor response to mean moisture in Field 2 is clearly anomalous at all 
sensor combinations with correlation ranging between 0.3 and 0,5 whereas 
correlation coefficients of the other fields are typically on the order of 
0.7 to 0.9. It should be noted that saturated conditions are included 
in the sample populations used to derive the correlation matrix given in 
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Figure 3.20. Backscattering coefficient at 4.6 GHz, 10°, HH polarization as a function of 
0-1 an gravimetric soil moisture. 




TABLE 3.15 


Correlation Matrix Between o°- „ and 0-1, 0-2 and 0-5 cm 

T,P )D 

Mean Gravimetric Soil Moisture for Each Field 


Angle Polarization Frequency 

(GHz) 


2.385 

3.291 

4.625 

5.625 

6.625 

7.625 

2.385 

3.291 

4.625 

5.625 

6.625 

7.625 



Linear Correlation Coefficient, R 


0-1 cm 

0-2 cm 

0-5 cm 

Field 

Field 

■ 1 f 1 - 1 — 

Field 

! i 



.898 .411 

.881 .406 

.849 .430 

.836 .514 


.805 .902 .547 .598 .808 .892 

.829 .892 .558 .672 .833 .894 

.876 .896 .390 .659 .871 .902 

.337 .877 .394 .534 .828 .873 

.881 .837 .418 .653 .879 .838 

.771 .825 .499 .732 .744 .828 

.744 .906 .532 .615 .746 .902 

.737 .786 .402 .606 .736 .778 

.740 .844 .546 .754 .741 .347 

.828 .805 .426 .595 .821 .797 

.743 .813 .444 .660 .722 .792 

.713 .699 .528 .722 .706 .704 






















Table 3.15 and that the linear correlation of o° with is slightly 
higher than shown for fields 3 and 4 when only nonsaturated moisture 
conditions are included in the analysis. Importantly, the exclusion of 
saturated data from field 2 does not significantly improve correlation 
coefficient. 

3.5.1 Poor sensor response to target inhomogeneity in Field 2 

The anomalously poor sensor response to mean moisture in field 2 
was examined with respect to the following potential error sources: 

1) failure of MAS system 

2) spatial mismatch of a° and M sampling 

3) temporal mismatch of radar and moisture sampling 

4) inadequate quantification of mean moisture due to non-uniform 

field distributions of: 

a) soil moisture 

b) bulk density 

c) soil texture 

Radar error was found to be inadequate to explain poor sensor response 
to moisture in field 2. External calibration of the MAS system with a 
Luneberg lens indicated no significant temporal drift of the system within 
the three month measurement period of the experiment (Dobson, 1979); and 
since all fields were examined throughout the entire period, all fields 
should be equally influenced by any calibration errors. Internal calibration 
of the MAS system was performed approximately twenty times within each 
data set by reference to a delay line. Inspection of delay line fluctuations 
indicated that the system was very stable within a data set; the delay 
line readings were typically found to vary by less than 1 dB. 

Analysis of the radar data with respect to the mean of sampled soil 
moistures assumes that the same spatially defined target is sampled 
contemporaneously by both the radar and "ground truth" efforts. Furthermore, 
it must be assumed that measurement errors associated with the radar and 
moisture sampling are random and normally distributed and thus would have 
minimal effect upon the sample means. During a radar data set of typically 
two hours duration the antenna boom was swept in an azimuthal arc of up to 
135° to cover the ground sampling grid shown in Figure 3.1. Since azimuth 


position was not recorded with the data stream, it was not possible to 
identify data sets where the entire field was not scanned. Partial scanning 
of a field would result in a sampling bias between the radar and moisture 
samples due to a spatial mismatch which would be data set specific. In a 
similar fashion, a sample bias between radar and "ground truth" data 
could occur when the two sets of target samples were not contemporaneous 
thus resulting in a temporal mismatch. 

Soil moisture is a time dependent variable. Moisture in the soil 
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surface layer was observed to change by as much as 0.12 g/cm in less than 
a 24-hour period during the 1975 Bare Soil Experiment. Temporal change 
in soil moisture can be particularly rapid when the soil surface is drying 
under the strong evaporative demands of solar irradiation and wind, 
especially when the soil is coarse textured or devoid of the moderating 
influence of a plant canopy. This temporal dependence of soil moisture 
poses a problem for the MAS system. 

A typical MAS 1-8 GHz data set requires 2 to 3 hours to collect an 
adequate sample size at all sensor combinations, while acquisition of the 
corresponding "ground truth" may require only 1/2 to 1 hour. Thus, the 
radar data includes a moisture related variation of 2 to 3 hours while the 
corresponding "ground truth" only accounts for a maximum of 1 hour 
temporal moisture variation. 

Three sets of procedures were incorporated in the 1977 data acquisition 
to help alleviate the uncertainty of correspondence between the soil moisture 
sampled by the radar and that sampled as "ground truth." First, the radar 
data acquisition time was reduced by limiting the number of incidence 
angles to three. Second, ground truth acquisition time was reduced such 
that soil moisture samples would reflect primarily spatial variations 
with minimal temporal variations. The use of additional manpower halved 
the acquisition time for ground truth; a normal ground truth data set did 
not include any sampling for bulk density, which can be a tedious and 
time-consuming process. And third, the temporal dependence of soil moisture 
was estimated by monitoring soil conductivity within the surface layers 
continuously through a data set. This was accomplished by using a set of 
conductivity probes buried in the soil. 
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Soil conductivity measurements were made at 15-minute intervals from 
three depth intervals at one location per field (Figure 3.21). The 
monitored depths were 0-1, 1-2 and 2-5 ck*. The resistance probes were 
coated with an abrasion-resistant porous dental plaster, Duroc, to a 
thickness related to the depth interval at which the probes would be inserted. 
The probes were emplaced in a given field shortly after field preparation 
and remained In equilibrium with the soil, unless vandalized as in the 
case of the probes at field 2, until the termination of the experiment. 

Figure 3.21 shows the probe placement; probes were oriented either 
horizontally or vertically within the soil in order to average moisture 
over the appropriate depth interval. 

Upon termination of the bare soil experiment, the conductivity probes 
were carefully removed from each field along with a block of the surrounding 
soil. The moisture dependence of the sensors was then calibrated in the 
lab. It should be noted that calibration curves are sensor and soil 
specific due to variance in the characteristics of the sensors themselves 
and also the variance between soil types in terms of hydraulic and salinity 
characteristics, Curtis and Trudgill (1974) offer an excellent review of 
the calibration technique, 

A representative calibration curve is shown in Figure 3.22. Although 
the laboratory calibration procedure was admittedly crude, it was still 
possible to define hysteresis effects on the sensors by the establishment 
of wetting and drying curves. The calibration curves for Field 3 were 
used to predict soil moisture from the mean conductivity recorded during 
"ground truth" acquisition. These moisture estimates compare very well 
with the mean measured "ground truth" moistures obtained for each data 
set as seen in Figure 3.23. The conductivity probes indicate changes in 
gravimetric soil moisture of as much as 5 l over three hour radar data sets 
at mid-day for drying soils. 

An example of short-term moisture flux can be seen in Figure 3.24. 

Two consecutive radar data sets were taken on Field 3; during this time 
the measured "ground truth" indicates that mean soil moisture in the 0-1 cm 
layer is fairly constant as a function of time but has an increasing spatial 
variance (especially toward drier conditions) while the 2-5 cm layer mean 
is increasing with an increasing spatial variance. The soil conductivity 
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Figure 3. 2 \ Placement of Conductivity Probes used to Measure 
Temporal Soil Moisture Response in the 1977 Bare Soil Experiment. 
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Figure 3.23 Conductivity Probe Estimate Accuracy In the 0-1 cm Layer 
Field 3 - All Data Sets Included. 


Drying or Wetting Phase is Determined by Rainfall History and the 0-1 cm Soil 
Moisture History. 




Figure 3. 24 Predicted Soil Moisture in the 0-1 and 2-5cm Layers for Two Data Sets 
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estimates of soil moisture, when plotted as a function of time, show that, 
at a given location, soil moisture is changing at a near constant rate. 

This change is not apparent from a comparison of the mean 0-1 cm "ground 
truth" moistures of the two data sets since the drying of the soil is 
masked by the spatial variance between the eight sampling locations. 

Due to constraints on time and resources, the other sensors were not 
calibrated with respect to the soils in which they were emplaced. The 
initial results however, indicate that this may be a very valuable 
technique to augment other forms of soil moisture sampling. It is highly 
recommended that future experiments include equipment to continuously 
monitor soil conductivity or resistivity; especially in view of the fact 
that resistance is proportional to soil matric potential and can be 
calibrated to gravimetric soil moisture by tensiometric field measurements 
and the laboratory determined desorption characteristics of a given soil. 

Since the conductivity sensors in field 2 were vandalized, it was 
not possible to compare its temporal variation in soil moisture within 
a data set with that of the other fields as estimated from the temporal 
change in soil conductivity. However, it is expected that this variance 
would be too small to account for all of the scatter observed in the 
radar response to field 2 (Figure 3.20). 

A much more plausible explanation of the anomalous behavior of field 2 
rests in the assumption that this field had a uniform distribution of 
moisture as sampled by both the radar and the ground sampling grid. If 
the true spatial distribution of target parameters such as texture, bulk 
density and moisture are not homogeneous but have large variance within 
the test field and complex spatial distributions, then the simple 
arithmetic mean moisture may be a very poor descriptor of the true mean 
moisture. The resulting bias imbedded in the mean would be compounded by 
any spatial mismatch between radar and "ground truth" for a given data set 
and also by any time transient conditions which could change the character 
of surface soil structure. During the course of the experiment., field 2 
was subjected to repeated partial flooding with water standing approximately 
6 cm deep over an irregular area of 30 to 50% of the test field. This 
produced a slaking of the soil surface in the flooded areas and resulted 


132 


in a puddled condition which increased soil density and altered the surface 
layer particle size distribution thus producing relatively impermeable 
surface crusts as a function of time. 

The surface layer distributions of soil texture and mean bulk density 
are plotted for each field in Figures 3.25 and 3.26 respectively. The 
spatial variance of soil texture as determined from the eight sampling 
locations is observed to be small for fields 1, 3 and 4 with gentle gradients 
(Figure 3.25). In contrast, field 2 displays high variance and rather 
steep textural gradients which divide the fields into regions composed of 
approximately 9 % sand and 35% clay separated by a region of 15% sand and 
29% cl ay . 

The mean bulk density distributions as plotted in Figure 3.26 for 
the 0-1 cm layer indicate that fields 1 and 3 are uniform compared to field 
2 which shows a distinct progression toward higher mean density for the 
portion of the field most subject to flooding and puddling effects. Field 
4 was not plotted because of the small sample size from which the mean could 
be determined. 

Prior measurements have indicated that both density and texture have 

profound effect on radar response to soil moisture due to the influence of 
these characteristics on the complex dielectric constant of moist soils. 

The within field distributions of mean soil moisture as computed from all 
data sets (Figure 3.27) exhibit spatial patterns which are dominated by 
variance in surface elevation and drainage characteristics (Figures 3.28 
to 3.31) and appear to be only weakly controlled by the spatial distributions 
of soil texture and soil bulk density. 

Further analysis of the effects of soil texture on radar backscatter 
response to soil moisture requires that the data input into the analysis 
be from homogeneous test fields with strong sensor response to gravimetric 
moisture. The erratic nature of the radar response to moisture in field 2 
is not consistent with prior experimental results and the data is therefore 
suspect. The observed non-uniformity of the surface layer of field 2 in 
terms of texture and bulk density combined with the compounding effects 
of sample bias due to spatial, and temporal mismatch between radar and 
"grot/nd truth" is believed to be responsible for the observed scatter in 
Figure 3.20. In light of the foregoing and the fact that fields 2 and 3 
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Figure 3.26. Contour plots of 0-1 cm Mean Bulk Density within Fielas 1, 2 and 3. 
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Figure 3. 29 Field 2, silty clay loam, surface elevation above an arbitrary datum plane. 
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Figure 3. 3i Field 4, silty clay, surface elevation above an arbitrary datum plane. 




are both classified as silty clay loams, the data from field 2 is not 
required for analysis of soil textural effects. 

4.0 SOU MOISTURE INDICATOR ANALYSIS 

Thus far, discussion of the data acquired during the 1977 Bare Soil. 
Experiment has facussedon radar response to gravimetric soil moisture. 

The linear correlation results presented in Table 3.15 show strong 
correlation between radar backscattering coefficient and gravimetric soil 
moisture for homogeneous fields of sandy loam, silty clay loam and silty 
clay (fields 1, 3 and 4 respectively). However, prior experimental results 
as described in chapter 2 clearly point out the failure of a gravimetric 
moisture indicator to account for the microwave sensor response to a 
complex scene of mixed soil textures. 

This chapter will discuss the 1977 Bare Soil Experiment results with 
the purpose of isolating those soil moisture indicators which are 
independent of soil texture in terms of sensor response, The combined 
textural data base of sandy loam, silty clay loam and silty clay fields 
will be used to establish the linear dependence of o° on gravimetric moisture, 
volumetric moisture, moisture above some transition moisture, normalized 
moisture and soil tension. 

In the experiments analyzed by Schmugge (1976a, b) and Batlivala and Ulaby 
(1977) and. described in Chapter 2, normalized moisture was defined as a 
percent of estimated field capacity (equations 2.39 and 2.31). In order 
to avoid the uncertainty of estimating moisture at a given tension as 
some function of other typically measured soil parameters (Table 2.10) 
for the 1977 Bare Soil Experiment, the soil moisture versus tension desorption 
characteristic curve was measured directly, albeit at two points. 

Triplicate core samples of the 0-5, 5-9 and 9-15 cm layers of each field 
were submitted to the National Soil Survey Laboratory in Lincoln, Nebraska 
for desorption of soil water with pressure extractors. The mean resultant 
values are presented in Table 4.1 for soil tensions of 1/3 bar and 15 
bars which are commonly used to describe the water retention of soil at 
field capacity and at the wilting point respectively. The 0-5 cm 1/3 
bar gravimetric water content of the undisturbed core samples was found 
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TABLE 4.1 


Soil Water Retention at 1/3 and 15 Bars Tension by 
Desorption, % Organic Carbon and Soil Bulk Density 
Determined by the National Soil Survey Laboratory, 
USDA Soil Conservation Service and Estimates of 
Water Retention at 31 Bars and Saturation. 










to vary from 14.6 for sanjly loam to 35.0 for silty clay. Also found 
in Table 4.1 are values of the water retention of air dry soil which 
approximates the hygroscopic coefficient of soil at 31 bars tension and 
the observed gravimetric water retention of the soil when a field was 
saturated and partially flooded. Table 4.2 demonstrates that the 
laboratory values of Table 4.1 are closely approximated by estimation 
algorithms from Table 2,11; the error of the estimate is observed to 
increase with depth below the surface, and Schmugge's approximation is 
generally found to be the more accurate of the two at 1/3 bar and the 
Dobson-Ulaby approximation more accurate at 15 bars tension. 

The measured values of moisture at 1/3 and 15 bars given in Table 4.1 
were used to define log-linear estimates of tension for moistures from dry 
to field saturation. These functions are presented and shown graphically 
in Figures 4.1 to 4.3 for the 0-5, 5-9 and 9-15 cm depth intervals respectively. 

A log-linear function was used to estimate soil tension instead of an 
exponential power function as proposed by Clapp and Hornberger (1978) because 
of 1) its simplicity, 2) a power function is not necessarily a good 
descriptor of the relationship between moisture and tension for soil textures 
finer than sandy loams, and 3) there would still persist error due to the 
hysteresis effect between adsorption and desorption curves which could 
not be accurately predicted (Mualem and Morel -Seytoux, 1978). In spite 
of its simplicity, the log-linear estimate still compares favorably with 
the desorption characteristic measured for various soils by Hoi tan (1968) 
and by Carlisle (1978). Figures 4.4 and 4.5 compare the log-linear 
estimate of tension for the 1977 Bare Soil Experiment sandy loam field 
with measured desorption curves of sands and sandy loams between .01 and 
15 bars. Figures 4.6 and 4.7 compare the 1977 Bare Soil Experiment silty 
clay loam and silty clay fields to measured desorption curves between .1 
and 15 bars. It should be noted that all soils plotted in Figures 4.4 to 
4.7 are from the surface horizon of cultivated fields with bulk densities 
and textures similar to the 1977 Bare Soil Experiment fields. 

Simple log-linear tension estimation algorithms, such as the ones 
shown in Figures 4.1 to 4.3, were used to compute estimates of soil tension 
from gravimetric moisture. A more complex model of the moisture-tension 
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TABLE 4.2 

Comparison of laboratory Determined Hater Retention at 1/3 and 15 Bars 
Tension with Empirical Estimation Algorithms Derived by Schmugge (1976a) 

and Dobson and Ulaby 


A = Schmugge's algorithms 

1/3 bar Mg = 25.1 - .21 (X sand) + .22 (X clay) 
15 bar Kg = 7.2 - .07 [% sand) + .24 (X clay) 


3 = Dobson and Ulaby algorithms 

1/3 bar Hg = 35.29 + 3.74 (X organic carbon) - .30 (X sand) - .15 (X silt) 
15 bar Mg = 4.8 + 2.47 (X organic carbon) + .24 (X clay) 
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Soil Tension (bars) 


Figure 4. 1 

1977 BARE SOIL 0-5cm WATER RETENTION CHARACTERISTICS 


•——Sandy Loam log (Tension) * 1. 8216-. 157447Mx 

• Silty Clay Loam log (Tension) * 3.9273-. 14634Mx 

* Silty Clay log (Tension) - 3.9071-. 125276MX 

Tension in bars 

Mx = 0-5cm Gravimetric Soil Moisture 
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Soil Tension (bars) 


Figure 4. 2 

1977 BARE SOIL 5-9cm WATER RETENTION CHARACTERISTICS 


• Sandy Loam log(Tension) ■ 1. 8216-. 157447MX 

A Silty Clay Loam log(Tension) - 4. 2023-. 165365Mx 

★ Silty Clay log(Tension) * 4. 441-: 141337MX 
Tension in bars 

Mx * 5-9cm Gravimetric soil moisture 



H6 


Soii Tension (bars) 


Figure 4.3 

1977 BARE SOIL 9~15cm WATER RETENTION CHARACTERISTICS 

• Sandy Loam log(Tension) = 1. 7187-. 129191Mx 

A Silty Clay Loam log(Tension) = 5. 7296-. 239659M* 

★ Silty Clay log(Tension) = 6. 2417-. 209322MX 
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Figure 4.4 

WATER RETENTION OF SANDS AND LOAMY SANDS 
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Figure 4.6 WATER RETENTION OF SILTY CLAY LOAM 
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characteristic assumes that all gravimetric moisture is desorbed at 10,000 
bars tension and inputs estimates of the water retention at 31 bars and 
saturation (0 bars) into the algorithm. 

For purposes of comparison, the prime soil moisture sensor configuration 
of 4 GHz, 10° incidence angle and HH polarization determined by Batlivala 
and Ulaby (1977) is used throughout most of the following analysis, 

4.1 Gravimetric and Volumetric Soil Moisture 

The 1977 Bare Soil Experiment radar response to gravimetric and 
volumetric moistures is shown in Figures 4.8 and 4,9 for fields of sandy 
loam, silty clay loam and silty clay. The least squares linear fits of 
o° with respect to moisture are shown for each soil texture with line 
lengths representing the range of measured moisture values. Individually, 
the pearson's product moment correlation coefficient between o D and moisture 
is quite high, ranging from =.7 to .9 at all depth intervals (Table 3.15*) 
and irrespective of soil moisture indicator. This emphasizes that for a 
given soil texture, regardless of the depth interval considered less than 
0-15 cm, radar backscatter at 4.6 GHz, 10° and HH polarization is highly 
dependent upon soil moisture. The strength of the dependence of a° on the 
soil moisture of a specific soil texture, as measured by correlation 
coefficient, compares favorably to the results obtained from the 1974 and 
1975 bare soil experiments (Batlivala and Ulaby, 1977). However, it is 
obvious from Figures 4.8 and 4.9 that the sensitivity of o° to gravimetric 
and volumetric soil moisture is highly dependent upon soil texture. 
Sensitivity is observed to be greatest for the coarse textured sandy loam' 
and least for the fine textured silty clay field; these are precisely the 
findings of Schmugge (1976a, b) and Batlivala and Ulaby (1977). When the 
three soil types are equally represented, the combined regression fits 
(Figures 4.8 and 4.9) have a low slope and poorly predict the observed 
variance between soil textures. The poor correspondence of the combined 
texture algorithm to the texture specific algorithms indicates that the 
correlation, of ct° with either gravimetric or volumetric moisture will be 
low for ground scenes of varied soil texture and will degrade rapidly with 
increasing depth interval. 
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Backscatter Coefficient o° (dB) Backscatter Coefficient o.° (dB) 


Figure 4.8 

Linear regression fits of radar response to aravi metric soil moisture 
in the 0-1, 0-2, 0-5 and 0-9 cm Layers. o° at 4.625 GHz, 10 , HH. 
1977 Bare Soil Experiment Data for fields of sandy loam, silty clay 
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4.2 Normalized Soil Moisture 


In equation 2.29, normalized moisture M p was defined as a percent of the 
1/3 bar moisture. However, the normalized moisture concept can be broadened 
to examine the effects of the tension at which moisture is normalized. 

M n = 100 x Mg/My (4.1) 

where: M p = normalized moisture, % 

My = Mg at tension, T 

Using the values from Tables 4.1 and 4.2 and the algorithms given in 
Figures 4.1 to 4.3, My can be varied as a function of T. At each My the linear 
dependence of a 0 on M p can be established and the resulting normalized 
moisture algorithms can be compared to those incorporating gravimetric, 
volumetric and Schmugge's percent of field capacity moisture indicators. 

The 1977 Bare Soil Experiment radar response to soil moisture expressed 
as a percent of the 1/3 bar moisture, which approximates field capacity for 
fine- textured soils, is shown in Figure 4.10 for the 0-1, 0-2, 0-5 and 0-9 cm 
depth intervals. Upon inspection, it is immediately evident that dependence 
of sensitivity on soil texture is greatly reduced when compared to gravimetric 
or volumetric moisture indicators. The remaining apparent dependence of 
sensitivity on soil texture was found to be statistically insignificant with 
the standard error of regression coefficients typically exceeding the 
difference between the sensitivities of various soil textures. Not surprisingly, 
the correlation between o° and % 1/3 bar moisture for the combined data base 
of all textures considered simultaneously is very high, in excess of 0.8 for 
the 0-1 cm layer, and degrades only slightly with increasing depth. 

A direct comparison of linear algorithms using gravimetric and 1/3 bar 
normalized moisture indicators is shown graphically in Figures 4.11 to 4.13 
for the 0-1, 0-2 and 0-5 cm depths respectively. As expected, the combined 
algorithms incorporating normalized moisture exhibit a dramatic improvement 
in correlation coefficient. Line lengths of the regression fits indicate 
the range of measured moisture values. In the gravimetric moisture plots 
{Figures 4.11a to 4.13a) the levels of n° for the 50, 100 and 150" of 
1/3 bar moistures are nearly constant for all soil textures, 
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Backsatter Coefficient o° (dB) Backscatter Coefficient o° (dB) 


Figure 4.10 

» 

Linear regression fits of radar response to % of 1/3 Bar AAoisture in 
the 0-1, 0-2, 0-5 and 0-9 cm Layers. o° at 4.625 GHz, 10°, HH. 
Data from 1977 Bare Soil Experiment for fields of sandy loam, silty 
clay loam and silty clay. 
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0-5cm % of 1/3 bar Moisture 

Comparison of 0-5 cm moisture algorithms for individual 
and combined soil textures at 4.6 GHz, 10°, HH polarization, 
a) gravimetric soil moisture and b) moisture normalized 
at 1/3 bar. 


For all three soil textures considered simultaneously and given 
equal weighting, the relationship between correlation coefficient and 
soil depth interval is shown in Figure 4.14 for moisture expressed 
gravimetrically, volumetrically, normalized to 1/3 bar moisture, 
normalized to 1 bar moisture and percent of field capacity as estimated 
by the hueristic relationship between texture and 1/3 bar moisture 
given by equations (2.30) and (2.31). The correlation with o° of both 
simple moisture indicators decreases rapidly with depth in spite of the 
fact that volumetric moisture shows some improvement over gravimetric 
moisture within the 0-2 cm interval, In sharp contrast, all of the 
normalized moisture indicators have correlation coefficients in excess 
of 0.8 for the 0-1 cm interval and correlations exhibit only mild decay 
as a function of increasing depth interval. 

The effect of normalizing moisture at the 1/3 bar tension on the 
0-N cm moisture sample distributions explain the improved correlation 
coefficient. Histograms of the radar response to the soil for each soil 
texture are shown in Figure 4.15 at 4.6 and 7.6 GHz at all angles of 

incidence measured during the experiment. When the cr° values of all 
textures are considered as one sample population, o° assumes a normal 
distribution at both frequencies and at all angles. The backscattering 
coefficient was also found to be normally distributed for the entire 
experiment at all other sensor combinations of frequency, polarization 
and angle. Histograms of the mean soil moisture sample populations for 
the 0-1, 0-5 and 0-15 depth intervals are shown in Figure 4,16. As soil 
depth interval increases gravimetric soil moisture assumes a multimodal 
distribution with the sampled mean moistures from each texture defining 
a statistically distinct sampling population, while for moisture expressed 
as a percent of 1/3 bar moisture the total sampling distribution becomes 
increasingly normal about a mean of approximately 95 % for the 0-15 cm 
layer. Comparing Figures 4.15 and 4.16, the a 0 distributions clearly 
display a high correspondence to the distribution of normalized moisture 
even on a field by field basis. 

It can be concluded from the foregoing that simple descriptors of 
soil water such as gravimetric and volumetric soil moisture, while adequate 
for single texture targets, are insufficient to describe the moisture 
perceived by radar for more complex and realistic targets. When soil 
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Figure 4.16 Histograms of 0-1, 0-5 and 0-15 cm Layer Gravimetric Moisture and Percent of 
1/3 bar Moisture for fields of Sandy Loam, Sitiy Clay Loam and Silty Clay. 

-- Sandy Loam 



0-lcm gravimetric moisture 0-lcm % of 1/3 bar moisture 
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moisture is normalized by the percent of field capacity concept, the 
observed soil moisture distribution becomes less dependent upon textural 
effects. Since the distributions of a° in Figure 4.15 are observed to 
be relatively independent of soil texture, the concept of soil moisture 
expressed as a percent of field capacity appears to be a good first 
approximation of the quantity of soil water to which radar is responsive. 
Furthermore, the normalized moisture indicators are superior to either 
gravimetric or volumetric moisture at all measured frequencies for 
ground scenes of mixed soil texture (Figure 4.17). 

My in equation 4.1 was calculated for each soil texture at a range 
of tensions between 0.04 and 31 bars using the simple log-linear func- 
tions from Figures 4.1 to 4,3. At each normalizing tension the following 
linear regression model was fit to each soil texture separately and to 
the combined data base: 

a° (dB) = a + b M r (4.2) 

The correlation coefficient of the combined texture algorithm as a 
function of the normalizing tension T is presented in Figure 4.18 for 
the 0-1, 0-2 and 0-5 cm depth intervals. It is evident that the optimal 
normalizing tensions, with respect to linear correlation coefficient, 
lie in the region between a: 0.3 and 2 bars of tension. Outside of this 
range the correlation coefficient is observed to decline, especially 
with increasing soil depth interval. It should be noted that the results 
presented in Figure 4.18 are based upon the full range of moisture con- 
ditions present during the experiment and includes data from fields 
which were 10 to 20% covered with standing water. 

While Figure 4.18 shows correlation coefficient to maximize at 
approximately 1 bar normalizing tension for all depth intervals (also 
Figure 4.14), the statistical significance of the apparent improvement 
in correlation coefficient at large depth intervals is not obvious. 

The predominance of high surface soil moistures prevalent throughout 
the experiment suggests that there should be only small contributions to 
a° from moisture at soil depths greater than 2 cm below the surface. It 
is argued that the high correlations for depth intervals greater than 
0-2 cm are the result of the high correlation between 0-N cm soil layers 
(Table 4.3). 
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Figure 4.17. Correlation coefficient as a function of frequency for 
various 0-5 cm soil moisture indicators at 10° incidence 
angle and HH polarization. Soil moisture is expressed 
as gravimetric, volumetric, % of 1/3 bar moisture, % of 
Schmugge's estimate of 1/3 bar moisture and % of i bar 
moisture. 
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TABLE 4.3 

Correlation Matrix of Soil Moisture 
Expressed as a Percent of 1/3 Bar 
Gravimetric Moisture for 0-N cm Depth 
Intervals, 1977 Bare Soil Experiment. 



Algorithm constants and linear correlation coefficients from equation 
(4,2) are plotted as a function of normalizing tension My in Figure 4.19 
for radar response at 4.6 GHz to the 0-1 cm moisture of all soil textures 
combined. The sensitivity of o° to M n of the combined texture algorithm 
(Figure 4.19) is compared to the sensitivities computed for each 
independent soil texture in Figure 4.20 as a function of normalizing 
tension. The under prediction of sensitivity by the combined algorithm 
at tensions in excess of 1.5 bars is largely responsible for the decline 
in correlation coefficient when moisture is normalized at tensions greater 
than 1 bar (Figures 4.18 and’ 4.19). The exclusion of flooded and * 
saturated soil data sets from the combined texture algorithm is observed 
to improve correlation coefficient by approximately 0,05 at all 
normalizing tensions (Figure 4.21). 

The empirical results obtained from the 1977 Bare Soil Experiment 
clearly demonstrate the superiority of normalized moisture over gravimetric 
and volumetric moisture as texture independent indicators of the soil 
water observed by radar between 1 and 8 GHz. However, the value of 
normalized moisture to a user community is inhibited by its poor 
interpretability over most of the range of normalizing tensions. Normalized 
moisture is perhaps most meaningful at the 1/3 bar tension where it 
approximates field capacity. When interpreted as a percent of field capacity, 
0 % is completely dry, 50% approximates the wilting point of agronomic 
plant species, 100% approximates field capacity of the soil and 150% is 
observed to correspond to complete soil saturation in the field, 

4.3 Transition Moisture Concept 

A transition moisture Wy was defined in Chapter 2 which separates 
the behavior of the dielectric constant of soil into two moisture response 
regions. The volume of soil water adsorbed to the soil particle surfaces 
which is structurally and dielectrically similar to ice is believed to 
be separated from bulk water by the transition moisture (Wang and Schmugge, 
1978). The tightly adsorbed water below the transition moisture should 
be several monomolecular layers thick depending upon the strength of 
the electric double layers. The transition moisture has been observed 
to vary with soil texture (Wang and Schmugge, 1978), and is found to be 
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0.01 0.05 0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 50.0 100.0 


Normalizing Tension (T) in bars (log Scale) 

Figure 4.20. 0-1 cm algorithm sensitivity as a function of normalizing tension for individual soil 

textures and all textures combined at 4.6 GHz, 10°, HH polarization. 
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from Fields 1, 3 and 4. 








smaller for sandy soils than for clayey soils as predicted by t>e double 
layer model and considerations of specific surface. 

The transition moisture concept was evaluated as a descriptor of 
the radar response to soil moisture by assuming that the bulk water in 
excess of the transition moisture is dielectrically equivalent regard- 
less of soil texture. The validity of this assumption should depend 
largely upon the salinities of the bulk soil solutions which the model • 
tacitly assumes to be equivalent. In addition, all measured M y during 
the bare soil experiment is assumed to be in excess of the transition 
moisture Wy. Thus, the postulated model of radar response to M y can be 
presented graphically as a series of parallel lines, one for each soil 
texture, offset from each other by their respective transition moistures 
and is described by the following equation for M y in excess of Wyi ■ 

o° = ctM y + 3-j F-j + + ^3 P 3 + **4^4 (4*3) 

where: a = sensitivity 

0 - intercept of ordinate axis 

F = dichotomous soil texture variable which can assume 
a value of 1 or 0 

1 = field 1, sandy loam 

2 - field 2, silty clay loam 

3 *• field 3, silty clay loam 

4 = field 4, silty clay 

Multiple linear regression fits of the above model were computed for 
M y as corrected for soil shrink-swell effects as described in Appendix A 
at 2.4, 4.6 and 7.6 GHz, HH polarization and incidence angles of 10° and 
20°. The multiple correlation coefficients R are between 0.7 and 0.8, 
even when erratic field 2 is included in the analysis (Table 4.4). R is 
observed to increase only slightly when field 2 is excluded from analysis. 
These correlations are comparable to those obtained from algorithms using 
normalized soil moisture indicators. 

The regression fits of equation 4.3, at 4.6 GHz, HH polarization are 
shown; in Figures 4.22 and 4.23 for 10° and 20° incidence angles respectively 
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TABLE 4.4 

Multiple Correlation Coefficients of Radar Response 
to Volumetric Moisture in Excess of the Transition Moisture 


Test Fields 
Incl uded 

Polarization 

i 

Frequency 

(GHz) 

1 

Incidence 

Angle 

.Soil D 
0-1 

1 

epth I 
0-2 

iterva 

0-5 

1 (cm) 
0-9 

All 

HH 

2.4 

10 

.745 

.722 

.752 

.702 




20 

.732 

.721 

.754 

.740 



4.6 

10 

.735 

.714 

.748 

.709 




20 

j 

.704 

.676 

.710 

.671 



7.6 

10 

.703 

.675 

.717 

.678 




20 

.679 

.647 

.689 

.658 

1, 3, 4 

HH 

2.4 

10 

.774 







20 

.734 





• 

4.6 

10 

.795 







20 

.744 






7.6 

10 

.746 







20 

.703 





Backscattering Coefficient, 
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Sandy Loam 

V/ 


Silty Cl; 
Loam 


Ity Clay 




Multiple Correlation Coefficient - .795 
Sensitivity « 4. 054 clB/0. 10g/cm 3 
Assumed 


0.1 0.2 0.3 0.4 0.5 0.6 0. 

0-1 cm Volumetric Moisture (g/cm 3 ) 

Multiple linear regression fit of 0-1 cm volumetric soil 
moisture in excess of the transition moisture with radar 
backscatter at 4.6 GHz, HH polarization and 10° incidence 
angle. 


Backscattering Coefficient , 


15 




-15 - Sensitivity ■ 3. 216 dB/0. 10g/cm 3 

Assumed 



O-lcm Volumetric Moisture (g/cm 3 ) 


Figure 4.23. Multiple linear regression fit of 0-1 cm volumetric soil 
moisture in excess of the transition moisture with radar 
backscatter at 4.6 GHz, HH polarization and 20° incidence 
angle. 





as a function of 0-1 cm volumetric moisture. The assumptions made in 

postulating the above model make impossible an exact determination of 

the transition moistures for each soil texture. However, the offsets 

of the regression fits along the abscissa indicate the relative differences 

between the transition moistures of each soil texture. These offsets, 

in terms of M v> are calculated by inverting equation (4.3) and solving 

for M v of each soil texture at a constant value of o°. According to 

the model, the relative differences between these values of M should 

be equivalent to the differences between VJ^ for each soil texture. By 

this method, W T of silty clay loam and silty clay is calculated to be from 

0.18 to 0.28 g/cm greater than of sandy loam depending upon sensor f, p and 0. 

Values of transition moisture can be converted to an estimate of 
the number of monomol ocular water layers enveloping the soil solids 
given the specific surface and bulk density of the soil. Specific 
surface was not directly measured for the soils examined during the 
1977 Bare Soil Experiment, but can be estimated from knowledge of the 
particle size distribution (Figures 3.2 and 3.3) and the clay mineralogy 
of the field soils. The volumetric moisture in one monolayer is 
approximated by: 

i 

M L * S x D w x p b x 10" 4 (4.4) 

where: M. = volumetric moisture of soil having one monomol ecular 

. L water layer, g/cm 3 

2 

S = specific surface of soil, m/g 

. 3 

P b = soil bulk density, g/cm 

D w = diameter of a water molecule 3 2.76A° 

Specific surface is estimated as the sum of mean specific surface values 
of the soil textural separates of sand, silt and clay given in Table 3.6. 

The specific surface of the sand and silt size factions are shown by 

o 

Figure 2.4 to be ”0.01 and 0.1 in /g respectively. The clay fraction of 
soils in the test site area commonly consist of approximately 50% 
montmoril Ionite and 30% illite (Dickey, 1979). The remaining 20% of 

the clay sized fraction is presumed to consist of non-clay minerals 

. 2 
with a specific surface of approximately 10.0 m /g. Net specific surface 
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is calculated using values of the specific surface of various clay 
minerals from Table 2.8 by: 

S “ + W fcm > + S cl a y^ c10y > U'V 

and 

S clay B °* 5 s m + °- 3 s i + °- 2 s 0 * (4-6) 

2 

where: S M » specific surface of montmorillonite =700 m /g 

Sj - specific surface of illite =100 m 2 /g 

2 

S Q « specific surface of non-clay minerals =10 m /g 

0 

Solving equation (4.5) for S using S c -j - 382 m /g from equation 
(4.6) yields the specific surface estimates presented in Table 4.5 for 
each soil texture (ignoring the specific surface of the organic soil 
fraction). Using the mean soil bulk densities in the 0-1 cm layer 
the volumetric moisture of one monomolecular water layer is com- 
puted from equation (4.4) to be 0,004, 0,048 to 0.046 and 0.063 g/cm 2 
for sandy loam, silty clay loam and silty clay respectively (Table 
4.5). The above values are very similar to the measured air-dry water 
contents of the soil (Table 4.1). Thus, the estimated values of 
are closely related to the hygroscopic coefficient at 31 bars of tension. 
If the offsets in volumetric moisture from the linear regression fits 
of o° to M v (Figures 4.22 and 4.23) are interpreted as the transition 
moistures Wy, then the ratio Wy/M^ yields a tightly adsorbed water layer 
on the order of 4 to 6 molecules thick which exhibits dielectric 
properties dissimilar to that of the bulk soil solution. This estimate 
compares favorably to the *8-9 molecular layers of water estimated to 
be below the transition moisture from dielectric measurements (Wang and 

Schmugge, 1978). Furthermore, the computation of estimated W- from the 
15 bar water retention of the soils examined during the 1977 Bare Soil 
Experiment by equation (2.13) yields values almost identical to those 
obtained from comparison of the multiple linear regression fits -in 
Figures 4.22 and 4.23. 

4. 4 Soil Tension 

The data acquired during the 1977 Bare Soil Experiment was used to 
test the hypothesis that radar response is linearly dependent upon soil 
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TABLE 4.5 

Estimates of Specific Surface and the Volumetric Moisture of One Monomolecular 
Water Layer for Each Field Examined During the 1977 Bare Soil Experiment. 


Field 

P b 3 

g/cm 

Percent 

Sand Silt Clay 

9 

Specific Surface, m / g 
Sand Silt Clay 

Total 

m l 

L 3 
g/cm 

1 

1.43 

55.0 

42.2 

2.8 

.0055 

.0422 

10.696 

10.7435 

.0042 

2 

1.22 

9.1 

53.9 

37.0 

.00091 

.0539 

141.34 

141 .39481 

.0476 

3 

1.20 

11.5 

52.2 

36.2 

.00115 

.0522 

138.284 

138.33735 

.0458 

4 

1.30 

9.0 

45.1 

.45.8 

.0009 

.0451 

174.956 

175.002 

.0628 


tension. Since soil tension was. not directly measured in the field duirng 
the acquisition of radar data, values of soil tension were estimated from 
the desorption characteristics of each field given in Table 4.1 and shown 

in Figures 4.1 to 4,3. The assumption that laboratory defined desorption 
characteristics adequately define the actual soil tension of structurally 
dynamic field soils is particularly poor for soils of high clay content 
such as silty clay loam and silty clay which commonly exhibit large hysteresi 
effects. Thus, the presented results can be considered as a "worst case" 
approximation which should show dramatic improvement if field tension 
values were available to account for the hysteretic dependence of tension 
on some function of moisture history. 

Using the simple log-linear estimates of soil tension from Figures 
4.1 to 4.3, o° at 4.6 GHz, 10° and HH polarization is plotted as a function 
of the logarithm of the 0-2 cm layer tension in Figure 4.24. The observed 
scatter in Figure 4.24 is felt to be primarily caused by hysteresis in the 
true moisture-tension characteristics of the fields which is not accounted 
for in the tension estimates used herein and caused by signal fading. 

While the 0-2 cm linear correlation coefficient r for the above sensor 
combination is -0.82 for data sets without partial flooding conditions as 
noted in Figures 3.8 to 3.11, the inclusion of all such data in the 
analysis only decreases r to -0.80. 

A more complex model of the estimated moisture-tension desorption 
characteristic is shown in Figure 4.25 for the 0-5 cm layer. The model 
assumes the following: 

1) all water is deso' hed at pressures of 10,000 bars (Brady, 1976); 

2) air dry soil moistures from Table 4.1 are equal to the hygroscopic 
coefficient at 31 bars (Lahav and Bresler, 1973); 

3) 15 and 1/3 bar moistures are those determined by the National 
Soil Survey Laboratory in Lincoln, Nebraska; 

4) saturation occurs at the values of Mg given in Appendix A as 
determined by the porosity of each soil, and 

5) tension is not allowed to fall below that at saturation which 
is defined at 0.001 bars for the very porous sandy loam and at ’0.1 bars 
for silty clay loam and silty clay which are fine textured and have small 

pore diameters. This is necessary because the. logarithm of zero , where ; 
saturation truly occurs , is infinite. 


c° = -2.849 log T (Tension 1+1.244 



0-2cm log (Tension), bars 

Figure 4.24. Radar response at 4.6 GHz, 10° and HH polarization to 
0-2 cm estimated soil tension as defined by a simple 
log-linear model of the moisture tension characteristic. 
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Soil Moisture Suction (bars) 



0 10 20 30 40 50 ' 60 


0-5cm Gravimetric Soil Moisture 

Figure 4.25. Desorption characteristic of 0-5 cm layer of sandy loam, 
silty clay loam and silty clay using estimates of water, 
retention at the ‘hygroscopic coefficient and at saturation 

I, fVvSF rS 


Linear regression fits of the radar data to the negative logarithms 
of tension as estimated using the above assumptions produce correlation 
coefficients comparable to those obtained using normalized moisture 
indicators. Figure 4.26 shows r to be approximately 0.8 arid constant as a 

function of depth at 4.6 GHz, 10° and HH polarization. Correlation 
coefficients at other sensor combinations are also similar to those for 
the normalized moisture indicators, 

5.0 CONCLUSIONS AND RECOMMENDATIONS 

Analysis of the 1977 Bare Soil Experiment substantiates an effect 
of soil texture on the behavior of soil water at microwave frequencies 
as noted by other investigators in Chapter 2. It has been shown that for 
realistic mixed texture ground scenes the simple soil moisture descriptors 
of total gravimetric or volumetric moisture fail to account for the inter- 
field variance observed by radar. Conversely, all moisture indicators 
which incorporate the concept of water potential have been shown to 
correlate highly with radar backscatter independent of soil texture. 

5.1 Evaluation and Interpretation of Soil Moisture Indicators 

The electric double layer model and considerations of the energy 
potential of soil water predict that as distance increases from soil 
particle surfaces there is an exponential decay of forces restricting 
the mobility of water and the capacity of water to oscillate at microwave 
frequencies . 

Evidence to date suggests substantial changes in the structure of 
water as a function of distance from hydrophilic soil particle surfaces. 
The density of water adsorbed to Ma-montmorillonite is known to be far 
greater than that of pure water up to distances of approximately 10 A 0 
(Figure 5.1) and its density is less than that of pure water, similar 
to ice, even at gravimetric moistures of 20-3055. Al though the transient 
structure of pure liquid water is only poorly understood, it is safe to 
assume that the high density of an absorbed water layer approximately 
3 to 4 molecular layers thick indicates a densely packed and highly 
ordered molecular structure. 

The tension of soil water is also known to decrease exponentially 
with distance from soil particle surfaces. Soil tension has proven 
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Correlation Coefficient (o° , Mx) 



Ficjure 4,26. Correlation coefficient as a function of depth 

interval for at 4.625 GHz, 10°, HH with various 
moisture indicators Mx for combined data from fields 
of sandy loam, silty clay loam and silty clay. 
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Density of Absorbed Water (g/cm 3 ) 



o DeWitt and Arens (1950) 
a Mackenzie (1958) 

• Mooney, Keenan, and Woods (1951) 
a Moorish (1954) 

■ Anderson and Low (1958) 


■P H ,o ^ 25° C 


Water Layer Approximately 3.0 A° Thick 


2. 0 3. 0 4. 

Water Content (gH 2 0/g clay) 

Figure 5.1. Density of adsorbed water on Na-montmorillonite (from Mitchell, 1976) 



itself to be highly interpretable and can be related to a host of sol 1 - 
water-plant relationships. Specific tension values have long been used 
in the literature to define critical points in the hydrodynamic behavior 
and agronomic influence of soil water (Figure 5,2V 

The volume of soil water held at a given total water potential is 
soil specific and is determined by the specific surface of a soil, the 
surface charge density, the concentration and type of soil solutes and 
the bulk structure and pore size distribution within the soil relative 
to a water table. Since surface charge density is known to be within 
the same order of magnitude for most soil mineral components, differences 
between soils in the volume of water held at high tensions as thin films 
enveloping soil particles is primarily determined by the specific surface 
of the soil. For most soils evidence indicates that the thicknesses of 
these water films, which include adsorbed and solvate water, are 
relatively independent of soil texture when measured in molecular layers 
of water (Figure 5.3). The volume of additional water which can be 
retained by soils at low tensions is primarily limited by the pore size 
distribution of the soil. 

Proceeding from the above and considering the soil water potential 
and dielectric constant responsive to an equivalent assemblage of forces 
existent within the soil medium, the following explanations can be 
offered for the failure of total gravimetric or volumetric moisture to 
account for the observed radar backscatter from diverse soil textures: 

1. Neither measure accounts for particle surface effects on 
adsorbed films of soil water which vary between soil 
textures and mineralogical compositions as a function of 
specific surface. 

2. Solute effects on the behavior of adsorbed and solvate water 
are disregarded by measures of gravimetric and volumetric 
water content. 

3. Variance between soils in the amount of water held at 
low tensions, which should be dielectrically equivalent 
to bulk water, are ignored by gravimetric moisture but 

i 

are considered by volumetric moisture because of the 
relationship between bulk density and porosity. Thus, 
volumetric moisture adequately explains the structural 
differences between various soils which control soil water 
retention at low tensions. 
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Hygroscopic 



Figure 5.2. 


Sketch of a typical soil -water systen. The wilting point and field capacity are 
defined at the tensions of 15 bars and 1/3 bar respectively (fron 1 . Wang and Schnugge, 



Figure 



Water Film Thickness (Number of Molecular Layers) 

5.3. Water potential as a function of the film thickness as 
measured by the weight of water retained divided by the 
weight of water needed to cover the soil with one molecular 
layer of water. The curves for clay, silt loam, and coarse 
sandy loam are brought close together. The silty clay loam 
is still quite different (from Taylor and Ashcroft, 1972). 
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Although solute concentration is known to affect the behavior of 
soil water even at high moisture values, up to M v * 0.5 for pure 
montmoril Ionite, it may be adequate in geographic regions of nearly 
constant soil salinity to use algorithms that predict the volumetric 
moisture in excess of a transition moisture. The moisture in excess 
of the transition moisture model as proposed by Wang and Schmugge 
(1978) is supported by the data from the 1977 Bare Soil Experiment. 

The transition moisture model was found to produce correlation coefficients 
from 0.7 to 0.8 for frequencies between 2.4 and 7.6 GHz at 10° and 20° 
incidence angles. 

The evidence to date suggests that the empirically derived transition 
moistures correspond to a water layer 4 to 9 molecules thick and soil 
matric potentials between 15 and 31 bars. Thus far, there is insufficient 
evidence from either soil physics or microwave studies of soil water to 
establish whether the transition moisture is defined by a sharp boundary 
differentiating two distinct phases of water each possessing unique 
structures and electromagnetic behavior or is defined at an arbitrary 
value within a transition zone. Studies of water density (Figure 5.1) 
suggest a sharp transition while electric double layer theory and studies 
of fluid mobility adjacent to soil particle surfaces (Kemper, I960; 

Kemper, et al., 1964; McBride, 1977) suggest a transition region. 

If further research demonstrates the transition moisture to be 
defined by the soil water retention at a unique tension such as 15 bars, 
then algorithms based upon this model will produce highly interpretable 
information to the soil moisture user community. 

Normalized moisture was found to be highly correlated to radar 
backscattar independent of soil texture. Using data from the 1977 
Bare Soil Experiment, linear correlation coefficients exceeded 0.8 when 
soil moisture was indexed to water retention at soil tensions between 
0.3 and 2.0 bars. While empirically determined algorithms using 
normalized moisture offer strong correlations, they may be heuristic 
in that the indexing procedure is not supported by a theoretical 
foundation. The estimating power of normalized moisture algorithms seems 
to reside in an inherent value of linking the quantity of soil water to 
soil tension T. Normalized moisture algorithms should be treated as 
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approximations of a more exact solution where for a given sensor 
configuration and surface roughness: 


'■/ 


T) 


(5.1) 


In spite of the observed predicting power of normalized moisture 
algorithms, they suffer from poor interpretability by the potential 
user community of soil moisture information. It will be difficult to 
define the physical properties of normalized moisture indices. For 
example, when moisture is normalized at a tension of 1/3 bar, the 
resultant index is only a crude approximation of the following: 

0% dryness 
50% « the wilting point 
100% « field capacity 
- saturation 


On the other hand, soil tension has long been used as the preferred 
descriptor of the soil-water system; its use has been inhibited by the 
difficulty of measuring soil tension under field conditions over a large 
range of values. Direct reading devices such as tensiometers are only 
accurate for wet conditions ranging from saturation to matric potentials 
of ^0, 85 bars, while non-direet devices such as soil resistivity blocks 
operate effectively only at higher tensions and must be calibrated for 
each specific soil sampled. Analysis of the 1977 Bare Soil Experiment 
shows a strong log-linear dependence of radar backscatter on estimated 
soil tension which is independent of soil texture. 

As a descriptor of remotely sensed soil moisture, soil tension 
seems to account for observed variance associated with the affinity of 
water for soil particle surfaces below the postulated transition moisture 
and account for variance due to soil structural control of porosity and 
the volumetric moisture content in excess of the transition moisture. 

The use of soil tension may also account for solute effects on the 
dielectric behavior of soil water; however, this has not been experimentally 
verified since only soil matric potential was measured during the 1977 
Bare Soil Experiment and not the total water potential (which would include 
the osmotic potential). 
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Furthermore, if the dielectric properties of moist soils can be 
treated solely as a function of soil tension, it can be argued that 
hysteresis in the moisture-tension characteristic will produce variance 
in radar backscatter which cannot be explained by either gravimetric or 
volumetric moisture even for a specific soil texture. Since the estimated 
matric potentials used in the analysis of the 1977 Bare Soil Experiment • 
are based solely upon laboratory desorption and not on in situ field 
measurements of tension (which would reflect hysteretic dependence, 
as seen in Figure 5,4), it is possible that hysteretic effects could 
explain much of the remaining scatter observed by radar in Figure 4.24. 
Until field measurements of tension are made in conjunction with micro- 
wave sensor observation of soils, the impact of hysteresis on sensor 
response cannot be fairly evaluated. 

However, if there is a hysteretic effect, then it should display 
itself most prominently in soils with high clay content. It is well 
documented that sandy soils exhibit very small hysteresis compared to 
clayey soils. In addition, soils that have a high percentage of expanding 
clay such as montmoril Ionite are particularly susceptible to hysteretic 
effects since episodes of soil shrinking and swelling act to modify the 
pore size distribution as a function of time. 

The above is supported by radar backscatter measurements from the 
1977 Bare Soil Experiment. Gravimetric soil moisture in field 1, a 
sandy loam, was found to be highly correlated with radar backscatter 
(in excess of 0,9) over the entire 3-month observation period; this is 
to be expected for a sandy soil with little hysteresis. Fields 2 and 3, 
silty clay loams, and field 4, silty clay, contained a high percentage 
of montmorillonite. In Table 3,15 the linear correlation coefficients 
between radar backscatter and gravimetric moisture typically ranged 
between 0.7 and 0.9 for fields 3 and 4 and between 0,4 and 0.6 for 
field 2; these comparatively lower correlations for high clay content 
soils are to be expected if there is a hysteretic dependence of the 
radar response to gravimetric moisture. The exceptionally low correlations 
for field 2 may be explained by the repeated flooding and dessication 
of the field which exacerbated changes in clay platelet orientation 
and porosity as a function of time which are typical of poorly drained 
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Moisture Content (cm 3 /cm 3 ) 



Suction (cm of Water) 

Figure 5.4* Hysteresis Scanning Curves. The experimental saturation 

boundary curve (a); the ultimate drying boundary curve (b): 
the wetting boundary curve (c); and a set of primary wet- 
ting curves (from Poulovassil is , 1973). 




montmoril lorn* tic soils. It is also significant to note that the radar 
backscatter response to field 2 was characterized by a pronounced time 
dependence explainable only in terms of soil surface layer structural 
changes and hysteresis. 

The foregoing discussion is relevant in that if there is a hysteretic 
effect on the radar response to gravimetric and volumetric soil moisture, 
then moisture estimation algorithms incorporating these soil moisture 
descriptors could be expected to exhibit poor performance even if the 
soil textures of a complex ground scene were mapped in detail and their 
locations and spatial extents were differentially treated in microwave 
sensor processing. Potential estimation error due to hysteresis is 
irrelevant for algorithms utilizing soil tension as the pertinent soil 
moisture quantifier. 

The major advantage of algorithms incorporating soil tension is 
that they express soil moisture with reference to the properties of soil. 
Total water potential accounts for the effects of soil particle surfaces, 
solute effects present in saline and alkaline soils of dry regions or 
the acid soil of bogs and highly fertilized agricultural soils, the 
effects of soil structure and hysteretic effects. Agriculturalists are 
interested in the effects of soil water on plants which respond to the 
energy status of water. Numerous studies of water uptake by plants 
have utilized matric potential as the relevant texture independent 
descriptor for the production of maximum yields of agronomically 
important crops (Table 5.1). Besides establishing the amount of soil 
water available to plants (Figure 5.5), soil tension determines (Table 
5.2): 

1. the amount of water that can be accepted by the soil 
before percolation starts, 

2. the aeration available to plants, 

3. saturation of the soil, 

4. the water available for drainage and run off, 

5. the vapor pressure of soil water if salinity is 
negligible (Figure 5.6) , 

6. the freezing point of soil water (Figure 5.7), 
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TABLE 5.1 

Matric potentials at which water should be applied for maximum yields 
of various crops grown in deep, well -drained soil that is fertilized 
and otherwise manages for maximum production. Where two values are 
given, the higher value is used when evaporative demand is high and 
the lower value when it is low; intermediate values are used when 
the atmospheric demand for evapotranspiration is intermediate. 

(These values are subject to revision as additional experimental 
data become available). (From Taylor and Ashcroft, 1972). 


Crop 

iMuiric potential 
{ joule f, kg) 

Eipiirolrnt 
matric vuctinn 
temibtirs) 

Reference 

Vegetative crops 
Alfalfa 

■ iso 

150 

S A Taylor' 

Beans (snap ami Inna) 

-75 to -200 

75 to 200 

Vittamet al. (1963) 

Cabbage 

-60 to -70 

60 to 70 

Vittani el al. (1963), and Pew 

Canning peas 

-30 to -50 

30 to 50 

(1958)’’ 

S. A. Taylor’ 

Celery 

-20 to -30 

20 to 30 

A. \V, Marsh,* and Marsh ( 1961 1 

Grass 

-30 to -100 

30 to 100 

Vissar (1959) 

Lettuce 

—40 to -60 

40 to 60 

A. W. Marsh.* Vissar (1959), and 

Tobacco 

-30 to -SO 

30 to 80 

Pew()958) b 
Jones et nl. (i960)’ 

Sugar cane 
Tensiometer 

— 15 to -50 

15 to 50 

Waterhouse et al. (1954)'* 

Blocks 

-100 to -200 

100 to 200 

Robinson (1963) 

Sweet corn 

-50 to -100 

50 to 100 

S. A. Taylor,* and Vittnm cl al 

Turfgrass 

-24 to -36 

24 to 36 

(1963) 

Morgan (1964) 

Root crops 
Onions 
Early growth 

-45 to -55 

45 to 55 

Pew (1958) h 

Bulbing time 

-55 to -65 

55 to 65 

Pew (1958)'’ 

Sugar beets 

—40 to —60 

40 to 60 

S. A, Tnylor* 

Potatoes 

-30 to -50 

30 to 50 

S. A Taylor.* Vittam et al ( 1963) 

Carrots 

-55 to -65 

55 to 65 

and Pew (1958) b 
Pew (195S) b 

Broccoli 

Early 

-45 to -55 

45 to 55 

Pew (195S) b 

After budding 

-60 to -70 

60 to 70 

Pew (I958) b 

Cauliflower 

-60 to -70 

60 to 70 

Pew (1958) b 

Fruit crops 
Lemons 

-40 

40 

A. W. Marsh* 

Oranges 

-20 to -100 

20 to 100 

Stolzyct al. (1963) 

Deciduous fruit 

— 50 to —80 

50 to 80 

A. W, Marsh,’ and Vissar (1959) 

Avocadoes 

-50 

50 

Richards cl al. (1962) 

Grapes 
Early season 

-40 to -50 

40 to 50 

A. W Marsh* 

During maturity 

<-100 

>100 

A. W. Marsh* 

Strawberries 

-20 to -30 

20 to 30 

A, W. Marsh,* and Marsh (196l i 

Cantaloupe 

-35 to -40 

35 to 40 

Marsh (1961), and Pew(l95S| b 
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Crop 

Matric potential 
(joules/kg) 

Equivalent 
matric auction 
(centibars) 

Reference 

Tomatoes 

-80to -150 

80 to 15U 

Vtttam ct al. (1958)/ and Viliam 

Daiian, is 

-30 to -150 

30 to 150 

el at. (1963) 
Schmeuli (1953/ 

Grain crops 
Corn 

Vegetative per km' 

-50 

50 

S. A. Taylor* 

During ripening 

-800 to -1200 

800 to 1200 

S. A. Taylor* 

Small grains 
Vegetative period 

-40 to -50 

40 to 50 

S. A. Taylor* 

During ripening 

-800 to -1200 

800 to 1200 

S. A. Taylor* 

Seed crops 
Alfalfa 

Prior to bloom 

-200 

200 

Taylor et al, (1959) 

During bloom 

-400 to -800 

400 to §00 

Taylor et si. {1959} 

During ripening 

-800 to -1500 

SQ0 to 1500 

Taylor et al.(1959) 

Carrots 

During seed year 
at 60 cm depth 

-400 to -600 

400 to 600 

Hawthorn (I95l)‘ 

Onions 

During seed year 
at 7 cm depth 

-400 to -600 

400 to 600 

Hawthorn (1951)* 

at 15 cm depth 

-150 

150 

Hawthorn (1951)* 

Lettuce 

During productive 
phase 

-300 

300 

Hawthorn ct al. (1956)’ 


Coffee Requires short periods of love 

potential to break bud dormancy, 
followed by high venter potential Alvin (1960) 


‘Unpublished research 

‘See also Marsh (1961) which includes n revision of few's work 

‘based on converting 50 percent available water to matric potential equivalent using Figure 14,16 
‘‘Tensiometers were located at IS inches, which is below the zone of maximum root concentration. 

'based on converting 50 percent available water to matric potential equivalent for loam to clay loam using 
Figure 14,16. 

'based on converting 67 percent available water to matric potential equivalent for clay soils using 
Figure 14.16. 

‘Resistance values of plaster moisture blocks were converted to matric potential from calibration of similar 
plaster units. 
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Soil Moisture Tension pF 

















TABLE 5.2 

Soil -Moisture Tension Relationships 
(from Kohnke, 1958). 
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pF Soil Water 


Figure 5 


Figure 5 
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0 10 20 30 40 50 60 70 80 90 100 
Relative Humidity 

.6. Relation between relative humidity and pF (from Kohnke, 1968) 



Freezing Point Depression, °C 

.7. Relation between the freezing point of water and pF 
(from Kohnke, 1968). 


7. the hydraulic conductivity of the soil, 

8. the evapotranspi ration rate during periods of high 
atmospheric demand (Figure 5.8), and 

9. when a soil is in the optimum moisture range for tillage 
With agricultural implements. 

5.2 Recommendations 

In light of the results of the 1977 Bare Soil Experiment the 
following recommendations are offered for further investigations: 

1. There exists a need for experimental verification of the 
relationship between the complex dielectric constant of 
moist soils and soil water potential. Such an experiment 
should strive to establish the dependence of the 
dielectric constant on soil tension, gravimetric and 
volumetric moisture for a wide variety of soil textures 
and clay mineralogies and seek a definition of the 
transition moisture in terms of a :pecific tension or 
range of tensions. The value of soil water potential 

in describing the dielectric behavior of saline soils 
should be investigated at low microwave frequencies. 

2. Radar experiments should be conducted for a wide range 
of soil textures where soil tension is measured in the 
field contemporaneously with radar data acquisition, This 
could be accomplished using a network of tensiometers 

to measure tension up to =0.85 bars and resistivity or 
conductivity probes to be calibrated to soil tensions 
at drier moistures. Such an experiment would serve to 
verify the stated conclusions of the 1977 Bare Soil 
Experiment and investigate the existence of hysteretic 
effects on the radar backscattering response to soil 
moisture. The investigation of hysteresis effects is not 
suited to the past format of limited time window airborne 
sensor experiments but should be pursued via long term 
experiments using ground based microwave sensor systems. 
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Figure 5.8. Relative evapotranspiration rate as a function of soil 

water potential for different potential evapotranspiration 
conditions. The E tp values for the curve are expressed 

in millimeters of water evapotranspi red in a 24-hour 
period (from Taylor and Ashcroft, 1972). 
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Appendix A: Physical Model to Account for Soil Shrink-Swell Characteristics 

The total pore space of the soil defines that fraction of the soil 
volume not occupied by soil particles. 

T « 100 x (1 - P b/» p ) (A.l ) 

where: T * percent pore space 

= soil bulk density, g/cm 3 
Pp = particle density 2.65 g/cm 3 . 

When the soil is saturated with water, the pore space is completely 
filled. Thus, 

M sat ' T (A ' 2 > 

where: M sa ^ s the percent volumetric soil moisture at saturation. 

3 

Assuming particle density equal to 2.65 g/cm , the mean bulk density 
values measured during the 1977 Bare Soil Experiment are used to compute 
the values of M sat presented in Table A.l. M gat is the total possible 
volumetric soil moisture that soil layers can attain unless the soil is 
permitted to swell in volume. Due to the prevalence of heavy rains before 
and during the experimental period, many of the measured soil moisture 
values were found to exceed Mg sat as given in Table A.l indicating that 
soil bulk density was dependent upon moisture content. 

Montmoril Ionite is the predominant clay mineral within the test site 
and is widely known for its capacity to swell. Test plots of silty clay 
loam and silty clay all exhibited extensive cracking of the soil surface 
as dessication caused soil shrinkage. Although the bulk density of these 
soils was measured in the field over a relatively narrow range of mois- 
tures near field capacity, density was found to be linearly dependent 
upon gravimetric soil moisture with a correlation of up to -0.74. 

Because of the dependence of density on moisture within fields of 
silty clay loam and silty clay, the calculation of volumetric soil mois- 
ture values from the mean measured soil bulk density was judged to be 
unacceptable. In spite of measured linear correlations of -0.6 to -0.7, 
the limited range and sample size of the measurements yielded a low 
significance to empirical moisture-density algorithms. In order to 
estimate volumetric soil moisture as a function of soil density, field 
and laboratory measurements of density were inserted into a simple 
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TABLE A.1 


Total Porosity and Gravimetric Moisture at Soil Saturation as Computed From Mean 

3 

Bulk Density and Assuming Particle Density Equals 2.65 g/cm . 

T = 1 00(1 - P b/ P p ) 

H 9sat “ T/P b 
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1.339 


1.353 

48.9 

36.2 

9-15 

1.592 

39.9 

25.1 

1.464 

44.7 

30.6 

1.437 


1.410 

46.8 

33.2 



physical model proposed by Fox (1964). Also of interest are the data 
and models presented in Berndt and Coughlan (1977), Jamison and Thompson 
(1967), Lawrence (1977) and Parker, Amos and Raster (1977). 

The model treats soil volume change as a three dimensional process 
at moistures below field capacity to account for the presence of crack 
systems and as a uni -dimensional process at moistures above field 
capacity where soil cracks are closed and further expansion can only 
occur in the vertical dimension. The following assumptions are made: 

1) Volume change is normal '.,id three dimensional from oven 
dryness until all of the vertical interstices of the soil 
are just filled with water and no further horizontal ex- 
pansion is possible. The swelling is similar to that 
observed by Berndt and Coughlan (1977) for semi and un- 
confined soil cores. 

2) Soil volume changes uni-dimensionally in the vertical 
dimension at soil moistures greater than those required 
to reduce the percent volume of entrapped air to 3%. 

This swelling is similar to that observed by Berndt and 
Coughlan (1977) for confined soil cores. 

3) The model is tied to values of soil moisture, bulk density 
and air capacity determined for undisturbed core samples 
from each field at a soil tension of 1/3 bar. 

4) The volume of entrapped air cannot be reduced below 3 %. 

The three-dimensional model is given by: 

P bx ’-Ps/(P S /A + M x p s + E l > 1/3 < A ' 3 > 

and the uni-dimensional model by: 

p bx = p s /(p s /A + M x p s + E 2> 
which reduces to: 

p bx = 1/? ^/ A + M x + E 2 ^ p s^ ( a - 5 ) 

3 

where: A = absolute density of soil = 2.65 g/cm 

3 

p bx = bu-| k dens i t y °f so1 ” 1 moisture g/cm 
p $ = bulk density at which all vertical interstices 
in the soil are assumed filled with water. This 
is assumed to occur at 1/3 bar tension. 

\ 

A- 3' 


M x * gravimetric moisture, g/g 

E-. « air capacity at 1/3 bar tension, cm 3 / cm 3 

3 3 

Eg ** air capacity at saturation * 0,03 cm /cm 

Figures A.l to A. 3 show bulk density as a function of gravimetric 
soil moisture for fields 2 to 4 respectively as computed from equations 
(A. 3 and A. 4). The three-dimensional model predicts density will 
decrease gradually with increasing moisture. When air capacity de- 
creases to 3 % all vertical interstices of the soil are assumed to have 
been filled with water and further swelling must occur uni-dimensionally. 

The figures show the uni-dimensional volume change to be quite dramatic 

3 

because each additional gram of water swells the soil by 1 cm , Figure 
A. 4 compares the volumetric moistures predicted by the model with those 
calculated for each field using a moisture constant mean bulk density. 

At low gravimetric moistures, M y from the model tends to be slightly 
greater than that computed from mean bulk density. At moistures 
approaching and in excess of saturation, the model predicts M y much 
lower than that computed from mean bulk density. 

Figures A. 5 to A. 8 show the volume percent of air, water and soil 
as a function of gravimetric moisture for fields 1 to 4 respectively. 
Field 1, a sandy loam, is assumed to have a constant bulk density as a 
function of moisture, this is firmly supported by the field bulk density 
measurements and the fact that this soil exhibited no surface cracking 
when the surface layer was dry. Fields 2 to 4 are each plotted using 
values derived from the physical soil shrink-swell model. 
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Bulk density as a function of gravimetric soil moisture 
for Field 2, silty clay loam, as predicted by the 
physical model. 
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Figure A. 2. Bulk density as a function of gravimetric soil moisture 
for Field 3, silty clay loam, as predicted by the 
physical model. 
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Figure A. 4. 0-5 cm volumetric soil moisture as a function of 

gravimetric soil moisture, comparison of results from 
the physical model to account for shrink/swell effects 
and results using a constant mean bulk density, (a) 
Field 2, (b) Field 3 and (c) Field 4. 
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Figure A. 5. Volume percent of air, water and stil as a function of 
gravimetric soil moisture for the 0-5 and 5-9 cm layers 
of Field 1, sandy loam. Values are computed with the 
assumption that mean measured bulk density is constant 
as a function of moisture content. 
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Figure A. 6. Volume percent of air, water and soil as a function of 
gravimetric soil moisture for the 0-5 and 5-9 cm layers 
of Field 2, silty clay loam. Values are computed from 
the physical volume change model to account for soil 
shrink and swell as a function of moisture content. 
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Figure A. 7, Volume percent of air, water and soil as a function of 
gravimetric soil moisture for the 0-5 and 5-9 cm layers 
of Field 3, silty clay loam. Values are computed from 
the physical volume change model to account for soil 
shrink and swell as a function of moisture content. 
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Figure A. 8. Volume percent of air, water and soil as a function of 
gravimetric soil moisture for the 0-5 and 5-9 cm layers 
of Field 4, silty clay. Values are computed from the 
physical volume change model to account for soil shrink 
and swell as a function of moisture content. 
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